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1 Introduction 

Quantum evolution of a few particles can effectively be computed by using the Schrodinger 
equation. Applying the same first principles to macroscopic systems with many degrees of 
freedom, however, is practically impossible. Macroscopic theories have to be developed which 
retain the relevant features of the original problems but are simple enough to be computa- 
tionally feasible. 

Due to thermal energy, lattice vibrations in metals can result in random deviations from 
the periodic background potential. This vibration is usually modelled by many independent 
harmonic oscillators. These are called phonons, and they can be considered bosonic particles. 

The interaction between the phonons and the electron is a collision process in which a 
phonon can be emitted or absorbed, subject to momentum and energy conservation. The 
number of phonons is not conserved, hence they are described by second quantization. 

Our goal is to show that the long time Schrodinger evolution of a quantum particle (elec- 
tron) in a phonon field can be described on large scales by the Boltzmann equation in c? > 3 
dimensions. We use two sets of variables; [x, t) stand for the microscopic space and time, and 
let 

(X,T) := {ex,et) 
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be the macroscopic variables. Here e is the scahng parameter separating microscopic and 
macroscopic scales and we will consider the £ ^ scaling limit. The velocity is unsealed. The 
quantum dynamics is given in microscopic variables. The resulting Boltzmann equation gives 
the evolution of the time dependent phase space density, Ft{X, V), in macroscopic variables, 
hence it contains only the macroscopic features of the evolution. We recall that the hnear 
Boltzmann equation with dispersion relation e{k) and collision kernel (riV, U) is 

drFriX, V) + Ve(V) ■ VxFt(X, V) 

= J [a{V,U)FT{X,U) -a{U,V)FT{X,V)]dU (1.1) 

= I a{V, U)Ft{X, U)dU - ao{V) Ft{X, V) (1.2) 

with the total cross section ao{V) := J a{U, V)dU. 

Acknowledgements. This work initially was a joint project with H.-T. Yau and several 
ideas presented here have been developed in collaboration with him. I would like to thank 

him for the invaluable discussions and encouragement through the entire work. 

Part of this project was completed during several visits at the Erwin Schrodinger Institute, 
Vienna, and at the Center of Theoretical Studies, Hsinchu, Taiwan. The author is grateful 
for the hospitality and financial support. 



1.1 Definitions 

For convenience, we fix a convention to avoid carrying factors of 27r along the Fourier trans- 
forms. We define the measure dx on R*^ to be the Lebesgue measure divided hy {2ny/\ With 
this convention the d-dimensional Fourier transform (usually denoted by hat) is defined as 

and its inverse 

fix) = J'-' fix) = [ fip)e^P-^dp . 

In particular, J^ae^^'^dx — Sip), where Sip) is understood with respect to the dp measure. 
This convention applies only to (/-dimensional integrals, in one dimension da still denotes the 
usual Lebesgue measure and /f^ e'^'do; = 2n6it). For functions ip ^ (^^(R"') we also define 
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the measure 6{ip{pj)dp supported on the level set {ip = 0} by 

F{p) ^ 
h'Ho) |V^(p)|" 

where da is the {d — l)-dimensional surface measure inherited from dp. 

The configuration space is a big box A = := [—^/2^TL/2,^/2^TL/2]^ C R'^ in c? > 3 
dimensions with Jj^dx = L'^. Its dual is A* = Af""' := {\^L~^ZY C R''. We use the shorter 
integral notation for the normalized sum on the dual lattice 

With this notation J^e'^'^'^dx — S{k) and Ji^,e^^'^dk — 5{x) where the delta functions are 
with respect to the measures dk and dx. In particular, 5{k) is the lattice delta function, i.e., 
5{k) — L"^ \i k = Q and is zero otherwise. 

Convention. Letters x,y,z will usually denote configuration variables in A C R*^; k,p,q, 
r, u, V, ^ stand for momentum variables in A* C R*^; s, t, a, /3, e, u, 77, r, 9, cu are real scalars and 
i,j, K, £, m, n are integers. The same convention applies to capital letters. 

For simplicity, we will omit the domains of integration in the notations. In general, the 
notation / is used for rf-dimensional integration with respect to the modified Lebesgue measure 
on R'^ if the L — > 00 limit has already been taken. Otherwise, the domain of integration is A 
or A*, depending on the integration variable. For one dimensional integrations the limits will 
always be indicated. 

We will take the thermodynamic limit L — > 00 before any other limit and our result will 
be uniform in L. The compact configuration space is a mere convenience in order to define 
certain expressions rigorously; for most purposes the reader can always think of R'' instead of 
A and A*. 

The electron is considered as a spinless nonrelativistic particle with state space Tie :— 
Lpg^(Ai). We denote the electron dispersion relation by e{k), k e R*^, i.e., the electron 
Hamiltonian, H^, is multiplication by e{k) in momentum space {k e A*). 

The state space of m independent phonons is V,^ := S ®Y=i ^IcA^l) , where S is the 
symmetrization operator. Their Hamiltonian, in momentum representation, is a multiplication 
operator by 

m 

H^hih, k2,..., km) = J2uj{kj) , kj e A* , 
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where uj{k) is the phonon dispersion law defined for k G R*^. 

The Fock space of the phonons is Tiph '■= ©^^g'^p/i' ^^'^ phonon Hamiltonian is 
Hph — (BmH^. The total Hilbert space, including the electron, is Htot '■— 'He <S) Ti-ph- The 
m-phonon sector consists of wavefunctions — '^'^{x; ki, /c2, ■ ■ ■ , km) ^V-e® "W^, where x 
is the electron coordinate. The elements of 'Htot are denoted by ^ = (^''")^^o- The total 
noninteracting Hamiltonian of the electron-phonon system is i^e <8) 1 + 1 (H) -f/p/i, which we shall 
write as + Hph for brevity. 

To define the interaction, we need to define the phonon creation and annilihation operators 
c\ and Ck as 



(c,.*)'"(x;/ci,/c2,...,/c„) := Vm + 1 *"*+^(a;; /c, /ci, /cg, . . . , fc. 



;i.3) 



cl^)"'{x-ki,k2,...,km) :-- 



J2 ^"^-'{x; ^1, , kJ6{kj - k) (1.4) 



(hat denotes omission). These operators satisfy the standard commutation relations, i.e., 

[ck,c\]=5{k-k') 



and any other commutator is zero. In terms of these operators Hp^ — J u!{k)Nkdk, where 
Nk :— clck is the number operator. 

The equilibrium state of the phonons is the Gibbs state 

-fph := Z-^ exp ( - (3Hph + J Nkdk) = Z'^ exp (^J [- Puj{k) + /^JiV^dfej (1.5) 
with inverse temperature (3 > and chemical potential /i. Here Z = Za is the normalization. 



I.e., 



^ — Tr^ph ( / [ ~ ^^^^^ ^ 



Nkdk 



such that Tr-^ ^^jph — 1. We define 



to be the expected number of phonons in the mode k. 
The interaction Hamiltonian is given by 



He-p := iA / Q{k) 

J A* 



dk, 



;i.6) 



(1.7) 



(1.8) 



as a multiplication operator on T^e, where Q{k) models the details of the electron-phonon 
interaction and A is the coupling constant. We will choose A = y/e and this is the weakest 
coupling that yields a nontrivial (non-free) macroscopic evolution on the time scale t = T/e. 

The full dynamics is given by the Schrodinger equation 

idt^t = H^t, (1.9) 

where 

is the full Hamiltonian. In particular, for Q{k) = \k\~^ we obtain the well-known Frohlich 
Hamiltonian describing the polaron. 

The Schrodinger equation can be written as an evolution equation for density matrices Fj 
on Htot- 

idtrt = [H,rt]. (1.10) 
For example, (|1.1Q ) and ( |1.9|) are equivalent if Ft is the projection operator onto the state 



The density matrix formalism is more general, and is appropriate to describe thermal states. 

We let 7e = 7e,05 Tr 7e = 1, be the initial electron density matrix; if we start from a pure 
electron state tpo, then 7e = 7e,o is the projection onto ipQ. The total initial density matrix is 

ro:=7e®7p/i- (1-11) 



The solution of ([TTOl) is Ft = e-'^^Toe'^" 



We are interested only in the evolution of the electron, i.e., we want to compute Tr-niTtO 
where O is an observable acting on He- We can first take partial trace Tr^.^^^ =: Trp^ to 
integrate out the phonon modes to obtain a density matrix on He- The phase space properties 
of such density matrix are described by the Wigner distribution. 

The Wigner distribution of any density matrix 7 on He is defined as 

W,{x,v) := j^e-'^-y^{x + ^-)dy , (1.12) 

or in momentum space 
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where hat on density matrices denotes Fourier transform in both variables. We use the conven- 
tion that the hat on functions defined on the phase space L^(A x A*) means Fourier transform 
in the first variable only. In particular, the Fourier transform of the Wigner distribution is 

W,{i, V) := ^ e~'^-^W,{x, v)dx = j(^v + - ^) . (1.13) 

For pure states, 7 = \ip){4'\, we have ■y{x,y) = 4'{x)4'{y) and W-y is equal to W^, the Wigner 
transform of the wavefunction ip, in accordance with | |EY2| |. 

We also define the rescaled Wigner distribution to describe large scale (macroscopic) prop- 
erties of the density matrix 

W^iX, V) := e-''W^{j, v) . (1.14) 

1.2 Assumptions 

We introduce the notation (x) := (x^ + 1)^^^. We use indexed constants Ci, C2, . . . to quantify 
estimates in our assumptions. 

For the electron and phonon dispersion relations, e{k), uj{k) > 0, we assume symmetry, 
e{k) = e(— /c), uj{k) = uj{—k) and 

\V^e{k)\<Ci{l + {kf-^) £ = 0,l,...,2c^, (1.15) 

\Vuj{k)\<C2{l + {kf'^) £ = 0,1,..., 2d. (1.16) 
We also assume that the functions k t— > $±(p, k) := e{k -|- p) ± uj{k) satisfy 

lim $±(p, fc) = 00 , (1.17) 

and 

< C3 < Hessfe$±(p, k) < C4 (1.18) 
(the upper bound follows from (|1.15|) and ( [1.16| )). By degree theory, these conditions imply 



that k I— > $±(p, k) has only one critical point. 
Let 

E±{V. 0, S) := {k : |$±(p, k)-~9\< S} (1.19) 
be a small neighborhood of the level set k \—>- $±(p, k) = 9. 
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We remark that conditions ( |1.15| )-( ]1.18| ) also imply the existence of two constants g,C > 
depending on Ci, . . . , C4 such that 

sup E±{p,e,6)nB{q,g) <C6g'^~^ (1.20) 

whenever 6, g < g. Here | ■ | denotes the Lebesgue measure of a set, and B{q, g) is the ball of 
radius g about q G R*^. We also need a certain transversality condition on two such sets. 

Transversality condition: There exist positive constants g, C5 such that for any 61,62, g < g, 
any pi,p2 € R*^ and any 6'i , G R 

C,6,62g''-^ 



sup 



E±{pi, 01, 61) n E±{p2, 02, 62) n B{q, g) 



< , , . 1.21 

\Pi -P2I 



In particular, elementary calculation shows that all these conditions are satisfied in d > 3 
if e{k) = \k'^ and HV^cuHoo is sufficiently small. This is our primary example. 

To ensure that the statistical operator is trace class, we always assume that 

m{uj{k) - > Ce > . (1.22) 

k 

Note that the number density function M is symmetric and it belongs to C^'^(R'^). 

We assume that Q{k) G C^'^(R'^) is real (since -ffg-p is self-adjoint) and sjTiimetric, i.e.. 



Q{k) = Q{-k) = Q{k) , (1.23) 
and it has a fast decay up to 2c?-derivatives 

^max jViQ(A:)| < C,{k)-''-^' . (1.24) 

The initial electron density matrix 7e = 7^ will depend on e so that it has a macroscopic 
profile in the scaling limit. We assume that the limit 

F^{X,V) ■= limiy4(X,y) (1.25) 

exists weakly in 5(R^'^). In addition, we assume that 

limsup / (p)^'^+^^7e(p,p)dp < Cg < 00 . (1.26) 

For example, 7g can be a pure state, 7^ := \ip'^){ip^\, with a WKB wavefunction 

^^(x) = e'^/2^(ex)e^^("^)/" , 

where A,Se 5(R'^). 



1.3 Main Theorem 



Theorem 1.1 Let A = ^/e and let solve the Schrddinger equation with initial con- 



dition Fq = 7g ®'-fph, where the initial electron density matrix satisfies l \1.2d^ ) and l \1.2(^ . We 



let := TrphTf be the electron density matrix at time t. 

We assume that the electron and phonon dispersion relations satisfy l \l.l^ - i\1.18{ ), ( \1.21\ ) 
and ( \1.22l ), while the interaction function Q{k) satisfies i . 2!^) and ( \1.24 )- Then for any 
T > 

lim lim W'e {X, V) = Ft{X, V) 



'T/e 



weakly in SilV^ x R''), and Ft satisfies the Boltzmann equation ^1.2J with initial condition Fq 
and collision kernel 



a{V,U): = 27r\Q{U -V)\''l^{N'{U -V) + l)6[e{V)-e{U)+uj{U -V)) 

+M{U - V)5(e{y) - e{U) - uiV -U))^. 
These two terms correspond to phonon emission and absorption, respectively. 



(1.27) 



Remark 1. The Bohzmann equation is irreversible, while the Schrodinger dynamics is 
reversible. This is not controversal, since the time evolved quantum state cannot be recon- 
structed from the solution of the Boltzmann equation. The macroscopic equation gives only 
a robust information on the evolution: it contains neither the phonons nor the microscopic 
details of the electron state. The former was neglected when taking partial trace, the latter 
is lost in the scaling (weak) limit. 

Remark 2. Scaling limit of the dynamics of a quantum particle in a weakly coupled random 

PLW|| and for long time in [[li)Y2|| . A similar 



potential was obtained for short time in 



result was obtained in ||EY1|| in a different limiting regime, in the so-called low density limit. 
Despite these similarities, the current proof differs substantially from ||EY2|| (see also a short 
announcement [0). We replace the idea of the partial time integration by a more effective 
classification of the indirect Feynman graphs. This enables us to expand the Duhamel formula 
further than in |pi)Y2| ] and still control the error terms (so-called indirect and recollision terms). 
In particular, for any n we will be able to gain an extra e*^ factor relative to the size of the main 
term with the exception of C{k)'^ pairings in the n-th order terms of the Duhamel expansion 
(see Propositions |7!^ , and |A.6|) . 



10 



Remark 3. The evolution problem for the electron density matrix 7^ = Tiph^t is formally 
equivalent to the Schrodinger equation 



with a time dependent Gaussian random potential V^{x,t) with covariance (in Fourier space) 

mM^)V^iq,s) = 5{p-q)\Q{p)f[[M^^^ (1.28) 

This means that the formal perturbation expansions of these two problems coincide term by 
term. However, this connection is only formal and in our rigorous proof we cannot and do not 
make use of it. 

Remark 4- Long time evolution of a microscopically localized electron weakly coupled to a 
phonon bath was studied in [|CEFM|] in the dipole approximation (see also references therein). 
The limiting equation is diffusive (Fokker- Planck) already on the first nontrivial time scale 
(in the van Hove limit). In this case there is no Boltzmann equation before diffusion emerges. 
The diffusion mechanism is quite different; it is a resonance effect between certain phonon 
modes and the eigenfrequencies of the confinement. 

Remark 5. In a more realistic transport model, electron-electron interactions should be 
included. This is a genuine many-body problem and is much more difficult. In classical 
dynamics an analogous result has been proven by O. Lanford |E| in the low density limit. 
There is no rigorous result in the quantum case. 

Convention. Throughout the paper we use the letter C to denote various constants that 
depend only on the dimension d and the constants Ci, C2 . . . quantifying the assumptions in 
Section |1.2| . For brevity, we usually neglect certain variables in the formulas if it does not 
cause confusion. In particular, we will omit the superscript e from Ff and 7^ . 



2 Structure of the proof 
2.1 Duhamel expansion 

We use the Duhamel formula (|2.1|) for H = Hq + He-p with Hq = He + Hph = e(— iV) + Hph. 
Recall that Hf. and Hph commute. 
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We are interested in the partial trace of Tt = e~**^roe**^ with respect to the phonon 
variables. In principle, we could fully expand e~**^ by Duhamel on both sides. But this 
expansion does not converge for large T = te, e := A^, unless we give a fairly detailed 
classification of Feynman graphs. For each graph of order n the best generally valid theoretical 
bound is {CX'^t)"'/n\ (although we essentially prove only (CaA^t)"/[n!]", < a < 1, in Lemma 
O), but then n\ is lost to combinatorics (number of graphs is ~ C^n\) and the series does not 



converge for t > C^^A^^, i.e., for T > C^^. Even if we gain an extra from most pairings 
(due to reabsorption and crossing), with some fixed k, it does not make the series converge. 

So the idea is to stop the Duhamel expansion once the expanded part is small enough. This 
means that either there are many absorptions-emissions or at least one reabsorption. However, 
immediate reabsorptions do contribute to the main term, and they have to be resummed (see 
later in ( |10.14| )). For simplicity, we use the word collision for absorption or emission. 

The Duhamel formula for H = Hq + H^-p states that for any Nq >1 

N=0 







/**[ds,]^°e-^^°^i/e-pe-*^i^°i/,_p . . . H,_pe-''^oH° . (2.1) 

•J 



For brevity, we introduced the following notation 

H* rt rt N N 



j=0 j=0 

where the star refers to the constraint t = Y,f=o ^j- Notice that the upper integration limits on 
the right hand side of ( p. 21) could be omitted because the lower limits and the delta function 
together guarantee them. 

For the threshold, we shall choose A^o = No{t) such that < A^^! < ^2.2 ^^^^ {^^)). 



Remark: This is not the choice in our earlier paper [ [li]Y2[| , where Nq\ ^ t, and we used 



partial time integration. Here we rather expand further to avoid partial time integration. This 
will require more effective control on classically irrelevant interference terms. 

Since we have to distinguish a few recoUisions, the Duhamel formula has to be written 



more carefully. Write i^e-p (|1-8| ) as follows 



He-p = ixj Q{k) [e-'^^-^cI - e'^'-^Ck] dk = iX J Q{k)e-'^-%kdk (2.3) 
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with 

6fc := 4 - c_fc, hence bl = Ck - c^^ = -&-fc , 

and we used the symmetry of Q. Here star means the adjoint, and bk, bl do not satify the 
canonical commutation relations, in fact 

] = foi' any k,m . 

We also introduce the notation 

n 1 

n Aj ■= A1A2 ...An and J] ■= AnAn-i . . . A2A1 (2.4) 

j=l j=n 

for the ordered product of noncommuting objects (operators) Aj. 

With these notations we rewrite ( p.l|) in the Fourier space of Tie as 

e-^*^ = E^^/jd.,]^ /(ndA:,)e--^»(ng(%)6.,e-^^^^ 

N=0 •' j=l j=l 

+A^° / [ds,]^' / ( n dfc,)e--^( n g(%)&.,e-^^o) . (2.5) 
2.2 Stopping rule for the Duhamel expansion 

The key observation is that the Duhamel formula (|2.5| ) is too rigid. There is no need to decide 
apriori how long we expand each term. Once a new step of the expansion is made, and a 
factor He_p added, it contains an integration over a new variable kj which is actually a big 
summation ( p.3|) . This means that each term ramifies into a sum of several new terms. We 
can decide individually for each term whether we want to continue the Duhamel expansion or 
stop. The stopping rule must depend only on the already expanded terms. We can stop the 
expansion for those terms which already had enough number of recollisions (it turns out that 
one is enough) or which already had many collisions. 

By recollision or reabsorption we mean kj + ki = for some j ^ i in (|2.5| ). It contains both 
physically relevant (e.g. cj.Cfc and c.fccl;.) and irrelevant (e.g., cj,c^^ and c^kCk) interaction 
terms. The contribution of the irrelevant terms will vanish in the L —>■ 00 limit. The recollision 
kj + fc^ = is called immediate if \j — i\ = 1, otherwise it is called genuine. The reason for 
this distinction is that immediate reabsorptions do occur without rendering the corresponding 
term small as e ^ 0. 
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Hence we will stop the Duhamel expansion if we either see A^'o collisions or if we see one 
genuine reabsorption. If we did not stop, it means fully expanded terms, which then has no 
reabsorption. 

We need several notations to define these terms rigorously. Let xi'^)^ v € A*, be the 
characteristic function at 0, i.e., x{v) = 1 if v = and x{v) = otherwise. 

Definition 2.1 Let Ai{n, N) be the set of (n + l)-tuple of nonnegative integer numbers m : = 
(mo, rrii, . . . , m„) with the property that n + 2\m\ = N, \m\ := J2j nij . Such an {n + \)-tuple is 
equivalent to an increasing subsequence fj, = /i(2), . . . , of the numbers {1, 2, ... , N} 

such that all differences //(j + 1) — //(j) are odd for all j — 0,1, . . . , N . For convenience we 
set ii{0) :— and ii{n +1) :— N + 1. The identification between m and /i is given by the 
relations + 1) — — 2m j + 1 for all j — 0,1, ... ,n. We will use both representations 
simultaneously. Let 

I,^I{rn):^{lJi{3) : j = 1, 2, . . . , n} = {j + 2 ^ : j = 1, 2, . . . , n} . 

i=0 

For each j — 0,1, . . . ,n let 

Ij := Ijim) := + 1, + 3, . . . , + 2mj - 1}, {with \Ij\ = ruj) , 

and 

I- ■■= I -{in) := {M , M + 2,..., ii{3) + 2m,}, {with \^\ = m,- + 1) , 
be its complement between two consecutive elements of the subsequence /i. Finally let 

n n 

J J{rn) := [j , r := [j , 

j=0 i=o 

and clearly {0, 1, ... , N} — J \J J'^ with J H J'^ — 0, \J\ — \rn\ and \J'^\ — n + \rn\ + 1. 

The increasing subsequence /x associated with ■m&A4{n, N) will determine the indices of 
those phonon interactions which are not parts of immediate reabsorptions, i.e., m encodes the 
immediate reabsorption pattern. Here N refers to the total number of phonon interactions 
and n of them are not part of immediate reabsorptions. In the sequel, these two numbers 
always have the same parity without further remark. The number mj denotes the number of 
immediate recoUision pairs between and + 1). 

14 



For any m G 7W(n, N) we define 

E{k,m):= n Xih + h+i) (2.6) 

66 J (m) 

to be the restriction onto momenta that respect the immediate reabsorption pattern given 
by m. Here k := {ki, . . . , /cat) stands for the collection of the phonon variables. The phonon 
momenta indexed by the subsequence n are called external, the rest are internal. The number 
of external momenta is n if 221 G M.{n, N). 
We then define the measure 

r.(m) 

d^ := / d^ 5(^, 2b) 



- dk:=j dk E{k,m) n (l - Xih + hO) , (2.7) 



where d^ := I\f=i dkj for simplicity. This measure excludes pairing among external momenta, 
while the internal momenta are consecutively paired. Notice that these restricted integrations 
make sense only for finite A. We also set 



/ 



#in,N) rim) 

dk:= / cl^' (2-8) 



m&Min,N) 

i.e., in this measure we exclude all recoUisions, which are not immediate. In the sequel, genuine 
recoUision will be simply called recoUision. 

Now we define the one recoUision terms. These terms will always be amputated, i.e., the 
last free propagator e~'*°^° will not be present. 

For any a e {2, 3, ... , n}, we let Main, N) C M{n, N) be the set of those m's such that 
mo = 0, fi{a) > 3 and if a = 2, we additionally require that mi > 1. The momenta with 
index ^(1) = 1 and /i(a) will form the (genuine) recoUision. We also let la = laini) '■— 
I{rn) \ {1; /^(c^)} be the set of the indices of the remaining external momenta. 

For m e Ma{n, N) we let 

r*{m,a) 

dk := / d^ x{ki + k^(^a))'^{k,rri) 

be the measure where the immediate reabsorption pattern given by m and a single genuine 
recoUision between the first and the /i(a)-th phonon is enforced. Let 



/*(m,,a) r , , 

dk:= j dkx{ki + k^^a))'^{k,m) [] {l-x{h + h')) (2.9) 



/ 



dk:=Y^ E / ■ (2-10) 

a=2 m&Main,N) 
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This measure expresses that there is no recoUision of external momenta apart from the desig- 
nated one between 1 and /^(a). 



Now we define the (amputated) operator acting on Titot 

N 

/O 



A{r,k,N) := A^/^"[d.,]f(nQ(^,)&.,e--^^") . (2.11) 



V%{t) := / dkAiT,k,N), V%{t) := E^n.ivM , (2-12) 



3=1 

We can integrate out this operator on different sets of k corresponding to various reabsorption 
patterns (superscript denotes the number of reabsorptions) : 

r#{n,N) N 

dkA{T,k,N), V%{t) := 

n=0 

Vl^ir) / dkA{T,k,N), V],{t) := E ^n,iv W ■ (2-13) 

•' n=2 

These are amputated objects. Depending on whether a term is fully expanded or not, we 
define non-amputated terms by letting e~'*^° or e~'^^^ act on it. Terms with a recoUision are 
always amputated. 
We define 

B{t,k,N): = A^/ [ds,]^e--°^°(nQ(fci)?^fe,e-''^''°) 

= i-i) r dso e-""^M(t - So, k, N) . (2.14) 
Jo 



The fully expanded terms are 

r* ...... , , _ r#i^,N) 

and we let 



ft f#in,N) 

^uMt) i-i) I e-'^'^^'V^it - So)dso ^ j dk B{t, k, N) , (2.15) 



N 



n=0 

The non-fuUy expanded terms are 



rt N 

KMi) ■■= H) / e-'^^^'KNit - ^o)dso , n%{t) := KMt) , 
Jo ' ' 



n=0 
N 



•'^ n=2 

We will use the following Duhamel expansion: 
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Lemma 2.2 For any fixed Nq > 1 we have 

No No-1 

e-^*^ = n%^{t) + Y: Kit) + ^Nit) . (2.17) 

N=3 N=0 

Proof. When expanding the Duhamel formula, at each step we generate a new fully 
expanded term and a term with one more interaction H^-p and a full propagator. Since -ffe-p 
contains a sum over all momenta ( p. 3D , we obtain |A*| new terms. Hence the expansion can 
be represented by a successively growing rooted tree-graph, where the vertices correspond to 
interaction terms from the expanded ife-p's. Each vertex ramifies into a terminal branch (free 
propagator) plus |A*| new branches whose endpoints are labelled by the momenta of the newly 
expanded interaction term. In this graph there is a unique path to the root from each vertex; 
the vertices on this path are called the predecessors of this vertex. In particular, every vertex 
(apart from the root) has an immediate predecessor, called its father. 

We say that a vertex forms a bond (immediate recollision) with its father if they have the 
same momentum label and if the father does not already form a bond with its own father. The 
notion of the bond is defined successively as the tree grows along further expansion. Vertices 
that are not part of a bond are called independent. 

We stop the expansion at any vertex if 

(i) it has A^o ~ 1 predecessor; or 

(ii) its momentum label coincides with the label of any of its independent predecessor 
different from its father (genuine recollision). 

These terms are included in and H]^, respectively, while the fully expanded terms 
(corresponding to terminal branches) are contained in the last term of (|2.17|) . □ 

Therefore we can write for any K > 1 integer 

7i = irp,,e i oe = [t) + 7^^ [t) , 

where ^ ^ 

lT''it)= E Trp.£^(t)ro[4(t)]* (2.18) 

N,N=Q 

corresponds to the main term containing no reabsorption and less than K collisions, and 

No-lNo-l No-lK-1 

iT--= E E Tr,,4ro[4]*+ E ETrp.4^o[4]* + Tr,,7^^,^o[no]* (2.19) 

17 



No r 1* / 1* r ^° 

+ H TiphHl^To T-C]^ + Trp/i S'^Tq + J2 '^n + '^No + X! ^at 



Af,7V=3 



Af=0 



N=3 



N=3 



2.3 Observables and Wigner transform 

Eventually we want to compute Tig 7^(9 with some electron observable O acting on Tig. For 
example, to determine the weak limit of W^^ in iS(R'^ x R*^), we have to test the macroscopic 
Wigner distribution ( |1.14| ) against a smooth function J{X, V) G i5(R°' x R°'). Let Je{x, v) : = 
e'^J{xe,v), hence Je{C,,v) := 6^'^J{^e^^,v), recalling that hat means Fourier transform in the 
first variable. 

A simple calculation shows that 

(J, M/;) := J J(X, V)W^^{X, V)dXdV = Tr^ 7^0, , 



where the observable is given by the kernel 



Oe{u,v) := JAu - V, 



U + V 



To control Trg '~ftOs we will need that O^O* is a uniformly bounded operator: 



sup II 0^0*11 < oo 



(2.20) 



(2.21) 



In our case it is enough if satisfies 



= sup / sup \Je{^,v)\d^ < oo (2.22) 

e J V 

since ||0e(9*|| < || J|P by an easy calculation. The estimate ( |2.22|) is satisfied for if 



/ sup |J(^,f)|d^ < oo 

J V 



(2.23) 



in particular if J G 5(R°' x R'^). 

In Sections we will show that for any T > 



lim sup lim sup lim sup lielK^ (^) 







(2.24) 



with t = Tie. In Section p!0| we will check that the weak limit of W^mainit-. satisfies the 
Boltzmann equation as L ^ oo, then e — > and finally K ^ oo. This will complete the proof 
of the Theorem. 
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3 Estimating the error terms: statement of the Main 
Lemma 



All error terms in Tif.'y'j^^ {t)0 have the structure 

TTe+phGit)To g{t) 



O 



with Q = S OT = H with appropriate indices (see ( p.l9| )). We use that Fq > and a simple 
Schwarz inequality of prototype 



Tt ATqB*0 < TiA TqA* + Ti 0*BToB*0 

1/2 T^*, 



.1/2 



= Tr ATqA* + Tr To'"B*00*BT^ 
< Tr ATqA* + Tr tI^^B*BTI^^ 

= Tr AToA* + \\00*\\ Tr BYnB* . 



Hence after eliminating the "cross-terms" in 7^''^(t) (see (|2.19| )) by Schwarz inequalities, 
we have the estimate for any 0<z^<l,-ft'<A^o 



( A'-l 7Vo-l 



+iy-'Non%(t)TJn%(t) 



N=K 



No 



iy-'NoJ2nUt)To[nUt) 

N=3 



(3.25) 



with Co ■■= C(l + sup^ II 0^0*11) < oo. 

The terms in the first two summations are fully expanded, their estimates are easier. 
The third and fourth terms still contain the full propagator. We estimate them in terms of 
amputated fully expanded terms by using unitarity and paying a price of an extra factor t. 
The prototype of such term is 

7(t) := Tiph dso dso e-^^°^A(t - So)ToB*{t - So)e''°'' , (3.26) 

where A, B denote one of the amputated ^ for 6 = or 6 = 1 and n < N < Nq (see 
( p.l2| )-( |2.13| )). By a simple Schwarz estimate inside the time integrations in ( |3.26| ) we obtain 



Tre7(t) <t-TT, 



dsoTiphAit - so)ToA*{t - So) + dSoTrphB{t - so)roB*{t - So) 



(3.27) 
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Therefore we can continue ( |3.25| ) using the definition of and a Schwarz inequahty: 



( K-l N No-1 N 

< CoTt^+J ^ E W E ^nAt)To[sl^{t)\ + E W E ^nAt)^o[sl^{t) 

L Ar=0 n=0 N=K ?i=0 



(3.28) 



n=0 







Ar=3 „=2 



The main lemma is the following: 

Lemma 3.1 Lei < a < 1 arwin < with the same parity. For simplicity, we let M := 
For the fully expanded objects: 



lim sup Tre+ph i"" jv(^)ro [i'°,jv(^) 

L— >oo \ 

f't/iis w// 6e Msec? /or N < K); 

lim sup Tr^+p/, i"" ^(t)ro K,7v(^) 

L— >oo \ 



2+\Af 



< 



p^ + (c^A^t)-|;,(iog*tr 



(3.29) 



< 



^ + (CA^t)-^(log*t)-- 



(3.30) 



(this will he used for K < N < Nq). For the amputated objects, N > 3, 



lim sup Tr,^pJv%it)To\v%it) 



L— >oo 



< - 



] ^^0f- + ] iCXHf^{log*tr'\iN > 7) ; 

(3.31) 



hmsup Tr,+pJvlj,it)To\Vlj,{t)r ] < ^ (CXHf 



nl 



(log* t)' + -{log* tr+''x{N> 7) 



t^ 



(3.32) 



We introduced the notation log* t := max{l,logt} and x is the characteristic function. 



Remark 1. Each estimate has two parts. The first terms include the contribution of the 
direct pairing term {"ladder" diagram) with possible immediate recollisions {"one-loop renor- 
malization") . The second terms contain the non-classical indirect terms (so-called "crossing" 
terms) . 

Remark 2. From the indirect terms in ( p.29| ) we gain only t^/^ instead of t, but the power 
of log* t is bounded. In ( 3.30| ) we gain full t-powers but we lose in the power of log* t. With 
more careful estimates it is possible to remove the logarithms and still gain full t-powers but 
we do not need it. 
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Remark 3. The basic idea in these estimates is that we do not gain factorials together with 
crossing or recoUision gains. This is only a technical convenience. In fact, instead of (|3.32|) , it 
is possible to prove 



lim sup Tie+ph Vl^^{t)ro te,jv(^) 

L—*oo \ 



*\1 [CaXHY n\ 



(3.33) 



t [M!]" t2 

which would be equally sufficient to estimate the recoUision terms. In this estimate we would 
gain a factorial together with a recoUision gain. We follow the first path since it is somewhat 
shorter. 



We use Lemma ^]T] to estimate the terms in ( p.28| ). We choose < a < 1 and 

2.2 logt 



hence 



log log t 



(3.34) 



for any 5 > 0, t > 1. 

First we let £ — then K ^ oo and finally z/ — > in ( |3.28| ). Recall that 
t = T/e, i.e., XH = T. Easy calculation shows that 



lim lim sup lim sup lim sup Tre7x (t) 

K^oo e^O L^oo 







e and 



(3.35) 



which proves (2.24) 



4 Formulas 



In this section is we derive expressions for the terms in Lemma 



4.1 Definition and estimates of basic functions 

We define the following functions for p,k E H^, a G {±}, rj ^ 0, a, s E H 

M{k,a): = \Q(k)\\^^ik) + ^) , (4.1) 

^aip,k): = e{k + p) + auj{k) , (4.2) 

e{s,p,n): = J2 [ e-''^'^-^P'^^+''^M{k,a)dk , (4.3) 
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We also define 



M*(A;) := max max \ViM(k,a) 

ct€{±} j=0,l,...,2d ' " 



and from (|1.24|) and tlie properties of JV from Section ^]2| we have 

\M*{k)\ < C{k)-^'^-^^ . 
Tlie following lemma estimates the functions G and T: 



Lemma 4.1 Under the conditions in Section \T73i and ( \4-(\ ) we have 

C 

snp\e{s,p,Q)\ < ■ 

sup |T^(a,p)| < C , sup \Trj{a,p)\ < 

p,a,r, V (a - e{p)) 



sup 



For ijj = {const.) we also have 



sup |V T^(a,p)| < C , sup |V T^(a,p)| < 
p,a,rj ^ V ^ (a - e{p)} 



< d . 



Moreover, the limit 



exists, and 



To+(a,p) := lim T^(a,p) 



/mTo+(a,p) = -7r ^ J \Q{k)\\Xik) + ^) 6{a - <^>^{p, k))dk . (4.12) 



ae{±} 



(4.5) 



(4.6) 



(4.7) 

(4.8) 
(4.9) 

(4.10) 
(4.11) 



Proof. Apart from the proof of ([4.12|) , we fix a and for simplicity we omit it from the 
notation, i.e., we assume that G and T are defined without the summation over a. The 
inequality ( [4.7| ) for |s| < 1 is trivial. For large |s| it follows from the stationary phase formula 



-is^{p,k 



'^M{k)dk\ < _(||(V)'^M|U. + ||(V)'^M|U. 



The prototypes of ( [4. 131) are the linear and quadratic cases; $(p, k) = u ■ k or $(p, k) = k'^. 
The linear phase factor is trivial: 

C 



-isu-k 



M{k)dk 



< 



{s\u\Y 



\{vrM\ 



(4.14) 
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For the quadratic case, by standard estimate on the Schrodinger evohition kernel 

C 



M{k)dk 



< 



(4.15) 



and ||M||ii < C|| (V)'^M||i2. Similar estimates are valid for other purely quadratic phase 
factors. 

The proof of ( [4. 13] ) for a general phase function uses a partition of unity separating the 
isolated critical point. In a small neighborhood of the nondegerate critical point, a smooth 
change of variables transforms the phase into a purely quadratic function and (|4.15| ) applies. 
Away from the critical point the gradient of $ is separated away from zero, i.e., (^4.14|) applies 
with |m| > C > 0. The size of the neighborhood and the Jacobians of the change of vari- 
able transformations are controlled uniformly in p using ( |1.18|) . The details of this standard 
technique are left to the reader (see also Sec. VIII. 2 of Pt|). 

For the first estimate in ( [4.8| ) we use the formula 

M{k)dk ^ 



a 



11 ds e* 

10 



3— is<E'(p, 



''^M{k)dk 



(4.16) 



$(p, k) + irj 

for any t] > 0. We estimate / e~^^'^''P'''^M{k)dk from ( |4.7| ) and use that d > 3 hence (s)"'^/^ is 
integrable. This inequality is a extension of Lemma 3.10 of [ [bL)Y2| . 

To obtain the factor decaying in {a — e{p)), we insert a smooth cutoff function k 0p{k) 
supported on the set {fc : |a — $(p, A;)| < 2}, so that 1— is supported on {/c : \a — ^{p,k)\> 
1}. Using (|TT3| ) and ( [TTBI ), we see that \V%{k)\ < C^{{pY + {kY) for any I. On the support 
of Op we repeat the argument above and use ([4.1 3|) with {'V)'^{M9p) < C{p)'^. On the 
support of 1 — the denominator can be estimated as 

C{k){{p) + {k) 



< 



c 



< 



\a — ^{p, k) + ir]\ {a — e{p + k) ± u (k)) {a — e{p)) 

again by (|1.15|) , and then the k integration is finite. 

We will sketch the estimate of the first term in ( [4.9| ), the other two are similar. We have 
from (p^) 



VpT^(a,p) 



< 



s e 



Vp<^{p,k) M{k)dk 



ds . 



Using stationary phase ( [4.13| ), the smoothness and decay properties of $ and M, we estimate 
the term in absolute value by C(s)^'^/^, which gives (|4.9| ) for c? > 3. 
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If C(j = (const.), then after a change of variables 



.^^M(i-rtd^ (4 17) 



Applying the same argument as in the proof of (|4.8| ) , using that V M is (i-times differentiable 
with decaying derivatives, we obtain ( [4.1(J| ). 

Finally, the limit ( [4 .111 ) exists by the estimate of the last term in ( [4.9| ). The formula ( [4.12[ ) 



is a straighforward calculation using the smoothness of e{k),uj{k),N'{k),Q{k). □ 



4.2 Graphical representation of collision histories 



We will write ^ ( |2.15| ) in momentum representation; pj stand for electron momenta, and k 



denote phonon momenta (see Section |4.4| later for more details). The history of the electron- 
phonon collisions will be represented by a graph, where the electron lines are solid and the 
phonon lines are dashed. The bullets represent collisions. We fix an orientation of the edges 
of the graph (Fig. 1). The momentum label on each edge expresses the momentum flow in 
the direction of the oriented edge. 

At each collision the adjacent three momenta are subject to momentum conservation. 
Fig. 1 represents a history with N collisions with initial momentum and final momentum 
Po- 



P P P P P 

1 2 N-1 N 



Observable I I A A A initial state 

H I ^2i ^ 3 I i^^N-l I 



Fig. 1. Electron (p) and phonon (k) momenta. 

r n * 

When we compute Trp/i£^° ^Fq ^ or similar quadratic expressions, we need another 
copy of the electron history and we will pair the 2A^ phonon momenta lines using the analogue 
of the Wick theorem. This will be represented by another copy of the graph on Fig. 1 but all the 
arrows reversed. The corresponding momenta will be distinguished by tilde. Pairing means 
joining the phonon lines and identifying the associated phonon momenta. More precisely, 
pairing of ka and fcb means identification ka = —k^ and similarly for two tilde variables. If ka 
and kb are paired, then the identification is ka = kb- The graph consisting of two copies of 
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Fig.l with a pairing of the phonon hnes will be called the pairing graph. The orientation is 
neglected on the figures. 

Immediate reabsorption occurs when neighboring phonon momenta are paired and these 
momenta will be integrated out separately. Therefore we will mainly focus on the pairing 
structure of the n external phonon lines. 

Definition 4.2 A pairing line in the graph is called an internal (immediate recollision) line 
if it pairs {ka,ka+i) or {ka,ka+i) for some a < N — 1. Otherwise it is called an external 
(genuine pairing) line. The skeleton of a graph is defined by removing all internal lines 
together with their vertices (Fig. 2). 



t • » 



I / 
I / 



y\ 



• — • — • — • — • — • — • — • — • — • — • — • — • — • 




Fig.2. A pairing graph and its skeleton (N=14, n=8). 



For most of our estimates only the skeleton will play a role and only the identity of the 
external lines should be kept in the notation. Hence we will relabel the indices of the phonon 
momenta (Fig. 3). Moreover, due to momentum conservation, the electron momenta between 
successive immediate recoUisions are the same so they can also be relabelled. 
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The size of the contribution of a specific pairing to TiphS^ j^Tq ^ is determined by the 
skeleton of the pairing. There is no recoUision hne, i.e., all external lines join an "upper" bullet 
with a "lower" one. The main (physical) contribution comes from the so-called direct pairing 
(or "ladder graph"). A pairing and the corresponding graph is called indirect or crossing 
if some {ka,kh) and {ka',kh') are paired with a < b, a' > b'. This notion is independent of 
relabelling. Fig. 4 shows a direct and an indirect graph. 



Fig.4. Skeleton of a direct pairing and part of an indirect pairing with a crossing. 

A graph is said to contain a genuine reabsorption if some ka is paired with kb, \a — b\ > 
2, or the same happens with some ka and k^. This concept is defined according to the 
original labelling, but notice that it depends only on one copy of the electron history. Fig. 5 
shows graphs with genuine reabsorption. Such graphs represent the trace of recollision terms. 



Notice that neighboring momenta may be paired in a skeleton if they have corresponded to 
a genuine recollision in the original graph and all the momenta in between formed immediate 
pairs (Fig.2). 



Fig. 5. Pairings with genuine recollisions (only part of the pairing lines are shown). 



4.3 The main representation formula 

Let n„ denote the set of permutations on {1,2,..., n}. 
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Proposition 4.3 For any n < N , N — n even, we have 



lim sup 



rr,+,j£:„%(t)rofe%(t) 



rn,m€M(n,N) 7i"6n„ 



(4.18) 



with 



Crn,m,^{i) ■= H (^'^T.iPn, Pn, t, k, a) Y^^^^^{t; p^, p„„ k, k, a) , (4.19) 

j = l,...,n 

where we define the measure 

n 

dUn{Pn,Pn,k,k,a) := ( n M{kj,aj)6{kj - k^Q))dkjdkj^dpndpJ{Pn - Pn)le{Pn,Pn) (4.20) 



on A* X A* X (A*)" x (A*)" x {±} c x R'^ x (R'^)" x (R"')" x {±}", and we let 



jr=o LfeG/; 



n 



n 0(sfe,Pi,^) 



fee/, 



(4.21) 



\^~si\i. n 



J|e^«'^[^fe)+^^] n0(5fe,P„fi, 



fee/, 



e'*'^ / da e-'*" \{ Rj'^' T^(a - fi,-, p,) 



(4.22) 



3=0 



where Ij,Ij etc. depend on m^m (see Definition |^. i|j and we define 

1 



:= Rj{a,Pj,VLj,ri) 
Rj := Rj{a,pj,nj,r]) 



a — e{pj) — Qj + if] 
1 



(4.23) 



a — e{pj) — fij — IT] 

In these formulas, pj 's are functions of pn and k = {ki, . . . , kn), and similarly for the tilde 
variables: 

n n 

~ ~ ~ (4.24) 



and 



Pj-=Pn+ Pj-=Pn+ 

e=j+i i=j+i 



= ^j{k,a) := Y ^i^i^e), Qj := Qj{k,ao n ^) = Y (^iT-^(i)^{ki>) (4-25) 
i=j+i e=j+i 



with a = (ai, . . . , a^). 
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Remark. This proposition shows that only those pairings are relevant that respect the 
no-recoUision rule. In particular, every pairing can be identified with a permutation on the 
(relabelled) indices of the external lines. In the sequel we use this identification freely. 

The following subsections contain the proof of Proposition ^^3. 



4.4 Duhamel formula in momentum space 



We express £^ (|2.15| ) in momentum space. Given k := (fci, ^2, • • • , k^) and pn, we define 



variables Po,Pi, ■ ■ ■ ,Pn-i as follows (Fig.l) 



N 

Pj=Pn+ E j = 0,l,...,iV-l. (4.26) 

i=j+i 



Sometimes these relations will be expressed as 

N-l 



f {Y[dp,)A{p,k) 



j=0 

with 

N-l N 

^ip,k) := n Kp^-P^ - ki), (4.27) 

j=0 i=j+l 

but mostly we consider Po,Pi, ■ ■ ■ ,Pn-i as functions of ki, . . . , k^ and p^- 

For expressions that are quadratic in E^ni ^.g. ( |4.18| ), we need another set of variables 



distinguished by tilde. Given k := (/ci, k2, ■ ■ ■ , k^) and piv, we define, for j = 0, 1, 2, . . . , — 1, 

N 

Pr-=PN+ E (4-28) 

e=j+i 



Note that the (fci, k2, ■ ■ ■ , kn), and {po,Pi, ■ ■ ■ ,Pn) variables in Proposition ^]3| will be only 
a subset of (fci, . . . , k^) and {po, . . . ,piy) after relabelling them (see later). 

Using this notation, we can rewrite the kernel B(t,k,N;pQ,pj^) of the operator B(t,k,N) 
( p.l4| ) in the Fourier space of He 

B{t,k,N-po,PN) : (4.29) 

„ N-l N 

= / ( n dp,)A{p,k) / [ds,]^e-'^«[^(^'°)+^-"l(nQ('^j)&fc,e-^^^[^(P^)+^^'^l) . 

This is an operator acting on Hph- Notice that po is not integrated out, i.e., there is a delta 
function S(^Po—PN~J2eLi kij on the right hand side showing that k,po,pN are not independent. 
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Equivalently, we can write 



N N 

B(t,k,N;pr, + Y.kj,PN) := / [dsj]^ e-''°^<P°^^''^'\]lQ{kj)bk^e-''^^<P^^^^ (4.30) 



if we consider pj (for j = 0, 1, . . . , — 1) as functions of /ci, . . . , and given by ( |4.26| ). 
We collect the terms depending on the phonon operators. Recall that 



-isHph J isHph _ -isuj{k) J 



-isui(k) J _ isu}{k) 



^-i{T-s)uj{k) _ ^i{T-s)u){k) 



Define 

then its adjoint is bl{s) = — 6-a:(s) and 

[bki'r),bmis)] = 5{k + m) 
Let 

G*{u,t) := {e-*™("W(M) + e*™(")(Ar(M) + 1)} , 

then 

^rphjphbu{T)b^{s) = TTph'yphbl{T)bl{s) = -G*{u, r - s)S{u + v) 
^^phlphbuir)bl{s) = Ti ph'jphbl{T)by{s) = G*{u, t - s)5{u - v) 

We define 

Tj := So + Si + . . . + Sj^i , fj- := So + 5i + . . . + Sj-i 

for j > 1 and tq = tq := and let r := (ri, . . . , tat), r := (fi, . . . , tat). 
The terms involving phonon operators in ([4.30|) are: 

/ N \ N 



-isoHp 



- i=i 



recalling the convention (|2.4|). We define 



Tik,k,T,f) : = Tiph 
= Trpft 



TV 



(4.31) 

(4.32) 
(4.33) 



(4.34) 



(4.35) 



l[bkA^i))e~''"'"'7phe''"'"^{ n b 

. j=l j=N 

IphblifN) ■ ■ ■ bl(f2)bl (ri)6fc,(ri)6fe,(r2) . . . 6fc^(r^) 



(4.36) 
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Combining ( p^.lSI ), ( [4.29| ), ( |4.35| ) and ( |4.36| ) we obtain an integral formula for the fully 
expanded term with no reabsorptions: 



(4.37) 



A 



X 



2N 



t* 



dkdk I [dsj]^[dsj]^T(fc, fc,r,r) / dp^dpN^ipN - Pn)%{pn,Pn 



N-l 



N 



j=o j=l 
( n dpj)A{p,k)e'''>'^P°^(^]le''^'^P^'^Q{kj 

j=0 j=l 

Here we used that Q{k) is real and that 6{po — Po) from Trg is equivalent to S{pn — Pn) 



4.5 Computing the phonon trace and the thermodynamic Umit 



We use Wick's theorem and the definition of ( f4.32 ) to compute the phonon trace ( [4 .361) 
appearing in (|4.37| ). We have to consider all pairings within the indices of the set {k,k} : = 
{ki, . . . , kiy,ki, . . . ,kiy}. We recall the definition of from ( p.8| ), and similarly we have 



dk 



rn£Min,N) 



(m) _ 

dk 



We also recall that m G Ai{n,N) is equivalent to a subsequence /i, which is viewed as a 
monotonic map /i : {1, . . . , ra} H-i> {1, . . . , A^}. The sets J, Ij, J depend on m and we use tilde 
for the sets / := /(m), Ij := /j(m) and J := Jim). 

Certain pairs are already prepared by imposing the immediate reabsorptions in the mea- 
sures via m^rn. We say that a Wick-pairing of {k,k} respects a given pair rn,rn& M.{n, N), 
if b is paired with 6+1 for all 6 G J and similarly for h E J (see ( ^■6[ )). 

If a Wick-pairing respects {m,m), then only the indices from I and / will be paired 
freely. Moreover, since the momenta with indices from / cannot be paired (see (|2.7|) ), and 
so are the momenta within J, the pairing occurs between the momenta kb {b G /) and k^ 
(6 G /). Hence any Wick-pairing that respects {m,m) can be uniquely identified with a map 
TT* G Vn '■= /in„yU~^ between the sets I and I. Here are the maps associated with m, m, 
and Vn = /iH„yU~^ is the set of maps of the form /i o tt o /i^^ : I I with tt G H„. 

The following lemma states that only those Wick-pairings are relevant that respect a given 
(m, m), the rest are negligible in the thermodynamic limit. 
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Lemma 4.4 Let m,m E A4{n,N) and let r, f be fixed times. For any function F of the 



momenta k,k and times r, f we have 



(m) Mm} _ _ 

dk I dk F{k, k, r, r)T(fc, k, r, f) 



Yl 5{kb- k^*^b)'jG*{kb,f^*(b) - n) 

b£l(rn) 



dk / d~kF{k,~k,r,f) 



'^(yk,T,m) X{k,f,rn) + O 



lAI 



where 

X{k,T,m):= n S{kb + kb+i)G*{h,n+i-n) (4.38) 

b£j{rn) 

expresses the trace of the immediate reabsorptions. Notice that the delta functions in X coin- 
cide with the characteristic functions ofB built into the measure /(^^ (13;. The error depends 
only on N and on the maximum of the function M . 

Proof. Consider a factor x{^b + ^fe+i) in the definition of H(fc, m). Suppose that the Wick 
pairing did not pair k}, with /cb+i. Then one can integrate out all the other /c, k variables in 
such order that eventually all delta functions disappear and we are left with an integration 

\H\\na 



dkbdkb+ixih + kb+i)H{kb, /cfe+i, r, r) 



< 



|A| 



(4.39) 



where H is the result of all the other integrations and clearly ||-ff||oo < C(^)||-P'||oo- 

Therefore only those Wick-pairings are relevant that respect the prepared immediate re- 
absorptions built into the measures via m, m and they are described by a map tt* G P„. For 
such pairing one can remove all further restrictions from the measure (imposed by (m), (m) 
superscripts in ( |2.7|) ) since S(fc, m) is imposed by the delta functions in X, while the comple- 
ment of the condition, that kb 7^ —kb> for b ^ b' E I (see ( p.7|) ), is small by a volume factor 
using an estimate similar to ( ^.39] ). □ 

Remark: Similar argument is valid for amputated and recollision terms, T>'^]^. Since we 
always take the thermodynamic limit L ^ 00 first, we can always assume that the prepared 
immediate reabsorptions are respected by the Wick pairings, modulo negligible errors. 



Using this Lemma and the fact that the expression in the big curly bracket in ( [4.37[ ) has 
a volume independent bound, we obtain from ( [4.37| ) 



Trg+p/i £^° jv('^)ro ^n,iv(^) 



E^. 

rn,m€M{n,N) TGlIn 



(4.40) 
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where for any m,mE Ai{n, N) and for any vr G n„ we define it* := fi o n o ^ ^ and 



<^bS{±} bei 



beiuJuj 



beJ 



J Q J Q - - 



b€lUJ 



(4.41) 



X 



N-l 



bei 

N 



n dp,)A{p,k)(l[e-''^<^^^ 

j=0 j=0 



be J 

beJ fcgj 

N-l ^ , , 

{l[dp,)A{p,m[e''^<^^^ 

j=0 j=0 



We used r^+i — Tb = Sb and the definition of M ( [4.1|) to combine the G* terms with the Q{k) 
factors: 

Q(r)2G#(r,C) = e'"'^"(")M(r, a) . (4.42) 

cr=±l 

4.6 Proof of the representation formula 

Recall that {0, 1, . . . , A^} is a disjoint union J U J'^. Let 6 t— > be the characteristic function 
of J, i.e., Xb = ^ b E J and Xfe = otherwise. Similarly, Xb is the characteristic function of 
J. 



Using the definition of tj, tj, we can rewrite the phonon phase factors in ( |4.41|) 



and 



with 



He 

bel 



b£j 



n ( n e 

j=o 6e/jU/= 



-istlflj{k,a_)+XbO-bi^(k 



■J~j' gio-bf^.(ft)a;(fcb) "Q ^10-555(^(^5) 



n( n 



6e/ 



gjs5[Qj(fc,CTo(7r*) 1)+X505w(fc5)]' 



fij := ^lj{k,a) := ^ abUj{kb) 

bel , b>ij.{j) 



Vlj := VLj{k, a o (tt 



(4.43) 



E %*)-i(6)^(^6) 

fee/ , fe>/"t(i) 

(recall that /i, /i are associated with m,rn). Here a := {ct;, : 6 G /} and usually we will omit 
the arguments. 
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Hence we have for m, l?l £ (^i N) that 



(4.44) 



be/uju J 



X 



beiuJ 



b&J 



N 



X [d.,]^ 



i=of,e/jU/= 



X / [d.~,]^ 



n n 

i=0fe6/jU/= 



In this formula we considered Pj,Pj as functions of pn = Pn and fc's or fc's, respectively (see 
( [4.26| )), and we used that 



Pb = P^i(j) + h for h e Ij 
Pb = Pt,{j) for be 1. 



c 

j ' 



and similarly for tildes. We also freely integrated out kb+i, ki_^^ for 6 G J, 6 G J since they do 
not appear any more in the formulas. 

Let p* := p* := p^Q), k* := A;* := a* := (t^q) for j = 1, . . . , n and p*Q := po, 

'■— Po- These will be the relabelled momenta (see Figs. 1-3.), temporarily distinguished by 
star. Then by (|4.26| ) we have p* = piy, = pN and 



P*j =Pn+ J2 P*j =Pn+ J2 K 

for j = 0, 1, . . . , n. Also notice that 



=i+i ^=i+i 



(4.45) 



e=j+i e=j+i 
Clearly Qj,Qj depend only on the star variables, i.e., with a slight abuse of notation we 
continue to denote them as 



Qj = Qj {k* , 01* ) , ^ j = ^ j {k ? QL* 



O TT 



(4.46) 



With these variables and after separating the terms with immediate recollisions, we obtain 
from (|4.44|) that 



2N ^ 

<T*S{±} 

j = l,...,n 



di^n{Pn,Pn,k*,k , a*) Y ~ Jt] p*^, p*^, k* , k ,a*) (4.47) 
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with 



-^teii} •' beJ i^j 
bejuj 



(4.48) 



X 



X 



[d^.]^=o n 



j=0 
n 



b€i'; be/, 



[dhZo n 



j=0 



ben 



beii 



where k* stands for the sequence kl, . . . and similarly for the other variables. We recall 
that p*, p*, rij, are functions of n, a*, k*, h and p* (m,m are considered given). 

There are two different expressions for ~ ^ after integrating out the internal momenta 
for all immediate recoUisions. Recalling the definition of (|4.3| ), we have 



t* 



[dsb]io n 

3=0 



n 



-iSble{Pj)+^j] 



(4.49) 



X 



[dss]f.o n 



3=0 



bel'! 



b&Ij 



Notice that all non-star variables have disappeared. 



The other expression for ~^ relies on performing all time integrals by using 



' m,m,TT 
N 



6=0 



da 



and we can regularize the ds integrations by shifting e{pj) ^{pj) ~ ^V- 
The result is 



t* 

[dsb]io n 



j=0 



n 

fee/"? 



n 



(4.50) 



/oo / 
dae-'*"n 
j=o LV 



nij+l 



n 



a - e{p*) - ilj +ir] J ^^^j a - e{p* + k^) - - abUj{h) + ir] _ 



and we will always choose t] := t ^. We can use a similar identity for the dsj integrals. 
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Again, we can integrate out kb, b G J, and similarly for the tilde variables. Recalling the 
definition of T ( [4 .41 ) we obtain from ( |4.49| ) 



j=0 



a - e{p*) - Qj+iri^ 



1 



nij+l 



(4.51) 



— e{p*) — — if] J 

and again, all non-star variables have disappeared. 

Finally, Proposition |4.3| follows from ( [4.40|) , (|4.47| ), ( [4.49|) and ([4.51|) if we remove the stars 
from each variable. □. 



5 Apriori bound 



We start with the following result that is weaker than (|3.29|) but we need it to estimate certain 
pairings later. 



Lemma 5.1 For any < a < 1, n < N , m,m E M.{n,N), tt G n„ tfe have 



l^m,m,7r(^)l ^ 



(recall that M := {n + N)/2), hence 



lim sup Tre+ph S^^^{t)To [i"" ^(t) 

L— >oo V 



2+\N 



< n! X 



(5.1) 



(5.2) 



Remark. The estimate ( ^.21) is summable over n and N only for short macroscopic time 
T = \H< Tn. 



Proof. The second inequality (|5.2|) easily follows from (|5.1| ) and ([4.18|) since the number 
of pairs m, m G (n, A^) is bounded by 

2 



n 



< C 



N 



(5.3) 



The bound (|5l|) is trivial if = 0. Otherwise we split \Y\ = |y |«/2|y |i-a/2 and 
we apply supremum bound in the first term. Using (|4.21|) and ([4.7|) , we see that 



feGJ 
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An easy calculation shows that {d > 3) 



t* 



1 +l-^"hi 



t 



M 



_ ^ 

(|J^|-1)!~ M! 



(recall that M := (n + N)/2 = n + |m|), hence 



(5.4) 



(5.5) 



For the \Y\^ part we use a Schwarz inequality to separate the tilde variables. The result 
from (|4:T9| ), (^ ), (g^) and (U) is 



< C'^'A^^'ItTt) sup / dz/:(p„,p„,fc,fc) 



M! 



(5.6) 



X • 



da Y\ \Rj 



j=0 



2-a 



3=0 



2-a 



with rj := t ^ and with a slight modification of the measure duT^: 

n n 

diy*{pn,Pn,k,k) ■= []\M*{kj)dkjdkj'j(^]\6{kj -k^(^j))^dpndpnS{pn-Pn)%{Pn,^ (5.7) 
i=i j=i 

and we recall the definition of M* from ( [4.5| ). The estimate of the two terms in the last line 
of ( ^.6|) are identical, so we consider only the first one. 

The case n = 0, mo > 1 is trivial by integrating out da, and collecting t™-oi2-a) _ ^m(2-a) 
using 



/oo 
-oo 



da 



<{Ct) 



mo 



This gives ( |5.1|) together with the t^^" = t—"- factor from ( ^.5|) . 

From now on we assume that n > 1. We can immediately integrate out ki, k2, . . . , kn and 
Pn- Moreover, we use the trivial estimate 



for all j. So we obtain 

<C^''A2^(^j t^^^^^l^'sup J df^{k,p^)y_Jal[\R, 

where 

dfi{k,pn) := 7e(Pn.,Pn)( n M* (A;j)dA;j) dp^ 



(5.8) 



j=0 



2-a 



(5.9) 



(5.10) 
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is a positive measure. The rest of the proof is similar to the proof of (3.16) in pi]Y2 
Fix all momentum variables, p, fc's, for a moment and let 

n-2 

Gia):=Y[\R,\. 

j=0 

For any given Pn,kn, Un we denote by q the probability measure on R 

da 



Z{pn, kn) a — A + ir] a — B + ir] 

with A := e{pn), B := e(p„_i) + cTnUj{kn) = e(p„ + A;„) + cTnUj{kn) and 

da 

-oo \a — A + ir]\\a — B + irjl 
(we neglect dependence in the notation). We have a simple estimate for Z: 

z< ^ 



re 

Z = Z{pn, kn) ■■= / 
J— ( 



1 + log+ 


A-B' 






V 



\A- B + ir]\ 
Rewrite the da integration in ( p.9[ ) as 

da Y[\Rj\=Z / ^(da)G( 

-oo • n J — OD 



3=0 



and use Holder (recall that 2 — a > 1) 

da n l^il)^ " < Z^'" r ^(da)G'(a)2-'^ 



— oo 



j=0 



/ 



/oo 
Q{da;pn, kn) 
-oo 

M*(A;„_2)d/c„_2 



(5.11) 



(5.12) 



(5.13) 



(5.14) 



We then obtain 

\Crn,^M < (^^')'^(m) ^^'""^'-'^^P / dpndkn[Z{pn,kn)t''%{Pn,Pn)M*{kn) (5.15) 

M*{kn-l)dkn-l 



2~a 



M*{ki)dki 



We have 



sup sup sup 

" 2. p„,kj+i,kj+2,---,k„ 



2-a 



X . . . X 



/ 



tt - e{po) + ir]- ELi (^ei^iki) 



2-a ■ 



M*{kj)dkj 



a - e{pj_i) + ir]- Y11=j c^e^ike) 



-;<Cat 



l-a 



2-a 



(5.16) 
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with T] := t ^. Recalling pj^i = pj + kj, the prototype of this inequality for 6 > 1 is 

f M*(k)dk ^ , 1 , , 

sup / - — ^ / — < Cbt^-^ . (5.17) 



p,0 



\6 - $±(p, k) + iri\ 

There will be several similar inequalities in the rest of the paper. They follow from the 
assumptions in Section |1.2| by elementary arguments. Their proofs use the same idea and the 
details will be omitted. Here we only sketch the proof of ( p. 17) ). 

Proof of (^■llj ). We consider a tiling of R'^ by identical cubes {Qj} of size g. We recall 



the quantity q from (|1.2CI|) . On each fixed cube Q = Qj we resolve the singularity on an 
exponential scale. We define the sets 

= Se{p) ■.= Qn{k : < |$±(p, k) - d\ < 2"^^} (5.18) 

for 1 < £ < 4 := [log* t], £ e N. We let 

Se, ■.= Qn{k : \<^±{p,k)-e\<2-''g} 

and ^0 ■=Q\iJi Si. The integrand on S^, £ = 0, 1, . . . , io is bounded by 2''^C(dist(g, 0))-^'^-^^ 
using the decay of M*{k). The volume of Si is of order 2^^ by ( |1.20| ). We can sum up these 
estimates for i to obtain Cqt^'^. Notice that for 6 = one obtains Cq log* t. Finally, the 
strong decay of the coefficients Cq ~ (dist(Q, O))"^*^"^^ allows to sum up the contributions 
from each cube. Notice that only the following property of M* was used in this proof 

Ell^1U-{Q.)<^- ° (5-19) 



We continue the estimate of (|5.15| ). Using (|5.16|) , we integrate out ki, k2, . . . , kn-i in this 



order, then we integrate out da using that g{da) is a probability measure: 



(5.20) 

We then use the estimate ( p.l3|) and integrate out dkn to collect one more factor Cat^~"' using 



( ^.1?| ). At the end, to do the dp„ integration, we use that / 7e(p,p)dp = Tr 7e = 1. The total 
power of if: is Ma + (2 — a)|m| + (1 — CL)n = n + 2|m| = using M = n + |m|. This completes 
the proof of Lemma |5.1| . □ 
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6 Fully expanded terms with small n and no recollision 



Here we prove ( |3.29| ). We use the representation ( |4.4(]| ). We can assume that n > 4, otherwise 



the apriori bound ( p.2| ) apphes. 

Definition 6.1 A pairing vr G n„ zs called crossing if tt ^ id. In particular, there exists 
a < b, with n{b) < n{a). In this case the pair of indices (a, b) is called crossing pair. 



Remark. This definition differs from the one given in | EY2|| (Definition 2.5). Here only 



pairing fines between tilde and non-tilde variables can form a crossing pair. 

The only non-crossing pairing is the so-called direct pairing, vr = id, it is estimated by 
]!]) and it gives the first term on the right hand side of ( p. 29 ). 



Lemma 6.2 Let n < N , m,m E A4{n, N) and let it G H.„ be a crossing pairing. Then 

limsupIC ~ < t-'/'{CXHf {log* tr. {Q.l] 



L—^QO 



From this lemma and the combinatorial bound (|5.3| ) the estimate (|3.29| ) follows immedi- 
ately. 

Proof of Lemma \6. 4 If (o, b) is a crossing pair, then for fixed b there could be many 



different a's which form a crossing pair with b. Consider the smallest one, i.e., 7r(6) < it (a), 
but for all c < a, 7r(c) < vr(6). This property will be called the minimality of the crossing pair. 
From ( ^191) , ( ^:22D , Q and the bound (U) we have 



/roo roc ^ 

dK{Pn,Pn,k,k) daY[\Rj\ dan 1^.1 • (6-2 



Here we considered Pj,Pj as functions of p„ = p„ and k or k, respectively, according to ( [4. 241) 
and recall the definition of Rj, Rj from ( |4.2I: 



We need the following estimate to take care of the a, a integrations. For any fixed index i 

(P^)<C^"(n(%))(Pn) (6.3) 

and the same estimate is true for (pi). This easily follows from {a + b) < C{a){b). 
We also define 

L{k) := {ky^M*{k) (6.4) 
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and we use that 

L{k) < C{k)-^''-^^ (6.5) 

by (p^. 

Let A := {a — l,a, and A := {7r(6) — 1, 7r(6), ra}. Notice that these are sets of three 
distinct elements. We then perform a Schwarz estimate in (|6.2|) to obtain 



dada 



(6.6) 



X 



X{\R,\X{\R,\=: {L) + {n.) 



where the decomposition (/.) + (//.) is according to the summation in the big square bracket. 

To estimate the term (/.), we use p„ and kj, j = 1, . . . ,n, as variables and we simply 
estimate the terms Ra, Ra-i by Ct each. Moreover, we can gain an extra factor (a) at the 
expense of {pnfYrj=i{kj)^ . Using (pj) and ([TTel) : 

Ct 

< ( a-jL^ — ■ " ^ ' ~" < ( r-f- 

(«) ,=1 



\Ra\ < 



{a - e{pa) - ^aik,a)) 
We also need the estimate 



{a) 



(6.7) 



(6.^ 



to insert a decay in ^^^{b) and we also insert an explicit extra IljX^i)"^ decay. Finally we 
estimate M*{k) by L{k){k)'^, which makes up for the Y[j{kj)^ factors used in these estimates. 
Then we can integrate out all kj and p„ freely, set p„ = p„ in Rn and we are left with 

f"^ dada 



(J.) < t-'mCXr''t'\^^snp f (]jLCk,)dk,) r ^ fd^^^^^^ 



x|-Rn| \Rn\ \R-K{b)-i\ \R-K(h)\\ n \Rj 



n -| 

n 



(6.9) 



Since Rj depends on kj+i, kj+2, ■ ■ ■ , kn and p„, the variables ki, k2, ■ ■ ■ , kn-i, kn can be 
integrated out successively in this order. The integrals of kj with j ^ 7r(6), 7r(6) + l successively 
eliminate the factors Rj-i that depend on pj_i = pj + kj, and they each give Ct by 

L{k)dk 



sup 



9 - $±(p, k) + it] 



< Ct 



(6.10) 
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(see ( p.l7| )) with k = kj, p = pj and 6 = a — Q.j, where we recall ( |4.46| ) that 



"J — X! ^n'^{m)^{kn 
m=j+l 



and that pj = Pn + J2m=j+i ^m- We recall that the function L{k) also satisfies ( p.l9D . 

The integrals of ^,^(6) and eliminate the factors and Rn(b) and they each give 

C log* t by the estimate 

L{k)dk 



sup 

p,0 



9 - $±(p, k) + irj 



< Clos*t 



(6.11) 



Moreover, from the A;7r(6) integral we gain an additional (a): 



L{k^ib))dk^(b) / C"(log*t)n,(A:,)' 



a - e{p^{b) + k^(b)) + ir/ ± uj{k„^b)) - fi^(f,) (p,r(f,))^ 
The prototype of this inequality is 



< 



{a) 



sup 



/ 



L{k)dk 



d - ^±{p,k) +17] (p) 



< 



C\og*t 



(6.12) 



with 6 = a — fln{b)- The proofs of (|6.11| ) and ( |6.12|) follow the same route as ( |5.17|) and we 
omit them. To gain (a), we use (|1.16| ) 



1 ^ c{n) ^ c-u,{k,y^ 



(6.13) 



{a - ^l) (a) {a) 

We insert these estimates into (|6.9| ); we have collected (Ct)"~^(Clog* t)^ so far and we also 
gained (a). 

Finally we integrate out a, a and p„: we obtain 

dada r le{Pn,Pn){Pn)^dPn ,g 

^ {a){a) J |a - e(p„) + i?7| |a - e(p„) + ir^l 

< {C\0g*tf j %iPn,Pn){Pnydpr, < {C log* tf 



using (|1.26|) . Altogether the first term in ( |6.6|) gives 

(I.) <t-^/^{CXHf (log* t)\ 



41 



We now estimate the second term (//.) in ( |S.6D after integrating out k/s and pn 

« n 

{II.) = t'/\CXYt'^-^^^Pj{UM*if'Mk^yPnleiPn,Pn) (6.15) 



X 



jGA -gj 



Recall that Rj depends on pj's that are functions of the variables kj and Pn- We express 
everything in terms of kj,pn variables by ( [4.24| ). Similarly, 

Pj = + km 

m : ■n(m)>j 

by the pairing kj = kT^^j). In particular we have 

Pa-l = Pa + ka = V + ka + kb with V := Pn + ' 

j>a + l 

Pnib)=u + ka with u:=pn+ ' 

j : 7r{i)>7r(!)) 

(l-irib) = (i* + (TaUJ^ka) with Cl* 1= ^ ajUj{kj) . 

i : 7r{j)>7r{6) 

Notice that f , m and Vt* are independent of fci, . . . , fca and kh- This is clear for f , and it follows 
for u and fi* from the minimality of (a, 6). 

• # ^ 



a - ^ / ' D 



P7r(b)-1 P7r(b) ="+ka 

Fig.6. Basic crossing estimate. 



Now we start estimating ( 6.15|) . First we estimate Rn{b)-i by Ct and we also gain a factor 
using 

= < Y[{kj)^ (6.16) 



a - e{p^b)-i) +irj- ^n{b)-i (") j=i 
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similarly to (|6.7D. We also need 



n(%)-'<c" 



analogously to ( |6.8| ). Hence we have 



(//.) < t'/\CXYt^\^\ snp f (f[L{k,)dk,) [ 



dada 
1^ 



(6.17) 



X / ^Pn{PnYle{Pn,Pn)\R^{b)\ \Rn\ \ H l^jf 11 l^il It-T^ H 77-^^ 



We then integrate out ki, . . . , ka-i using ( |6.10D and collecting This eliminates 

Rq, Ri, . . . , Ra-2- Notice that Rn{b) is independent of these variables. 
Next we integrate out ka- This occurs only in Ra-i and Rn(b) 

(Clog*t)2 



\Ra-l\ \R7r{b)\{Pa) L{ka)dka < 

Here we used the inequality 
sup 



L{k)dk 



\Pa -U\{a- Via) ' 

{C\og*tf 



6 - $±(p, k) + iTj 



< 



(6.18) 



(6.19) 



e -^±{u,k) + i7] {p)^ \p-u\{e) 

with k = ka, p = Pa, d = a — Qa, = a — Cl*. 

Proof of ( \6.1[^ ). We use the transversality condition ( |1.21|) and a similar resolution of 
singularities as in the proof of (|5.17 ). Again, we consider a tiling by cubes of size ~ g. We 
fix a cube Q and we recall the definition of the sets Si^p) ( [5.18| ). On the set Si{p) fl S-^u) the 
integrand is bounded by 2^~^^C{dist(Q,0))~^'^~^'^{p)~'^, the volume is bounded by 2^^^^^^C\p — 
u\^^. To ensure the decay in 6* for A; G Se{p), i > 1, we use (p)'^ < {dist{Q,0))'^{9)^^ by 
( |1.15| ) and (|1.16|) . On So{p) we have 



< 



(dist(Q,0))^ 



6 - $±(p, k) + irj {0) 

After summing these bounds for i,£ = 0, . . . , [log*t], then summing up the result for all Q, 
we obtain (|6.19|) . □ 



Continuing (|6.18| ), from {6) ^ = {a — Qa) ^ we also gain a factor (a) ^ similarly to ( |6.13D 
using the decaying factor Yijikj)'"^ in ( |6.17| ). Hence we have 



{II.) < ei\CtY-\C\og*t)\C\Ye\^\snp I ( n L{kj)dkj 



(6.20) 



°° dada 



n-1 



dPn{Pn)'^%iPn, Pn) \Ra\ \Rn\ \Rn\ H l^i 



.j=a+l 



\Pa - U\ 
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Let c > a + 1 be the smallest index such that 7r(c) < 7r(6) (see Fig. 6). Such index exists, 
if not else then c = b. By the definition of v, u, we have 



Pa — u = V + kb — u =: kc 



w 



where w depends only on fcc+i, ■ ■ ■ ,kn,pn by the minimality of c. Simply k^ is the momentum 
with the smallest index which appears in the difference 



V — u 



j>a + l j : 7r(j)>7r(6) 



i.e., the first momentum which appears in v but not in u; or if the smallest such index is bigger 
than b, then c := b. 

Now we can integrate out ka+i in ( |6.20|) to eliminate |-Ra|- This gives only a Clog* t factor 
by ( p. Ill ) if c > a + 1, and by 

L{ka+i)dka+i 



I 



a - e{pa+i + ka+i) +ir]± uj{ka+i) - fi^+i 
if c = a + 1. The prototype of this inequality is 

L{k)dk 



sup 

p,w,8 



— ^±{p, k) + if] \k + w 



kn+i +w\ 



< C\og*t 



< C\oK*t 



(6.21) 



and it is proven similarly as ( |5.17| ) . Here the point singularity around w also has to be resolved 
on an exponential scale; we find that within each cube Q the integrand is bounded by 2^"*"^ 
except for a set of measure ~ 2^^^^^^^^^ using ( |1.2CI| ). The details are omitted. 



If c = a + 1, then all the remaining factors |-RjP, j = a + 1, a + 2, . . . , n — 1, can be integrated 
out successively. The order of the integration is dka+2, • • • ? dA;„. This gives (Ct)"^*^^^, which 
makes altogether a t-power — 1/2, taking the t^/'^(Ct)^~^0—^ prefactor into account and 
recalling that n + 2\m\ = N. 

If c > a + 1, then we integrate out ka+2 ■ ■ ■ , ^c-i in this order (it is void if c = a + 2). 
These integrations eliminate |-RjP, j = a + 1, . . . , c — 2, they do not affect the denominator 



\Pa — u\ = \kc + w\, and they give (Ct) 
eliminate Rc-i- 



^ ^ hj ( p.lOj ). Now we perform the kc integral and 

L{kc)dkc 



a - e{pc + kc) +ir]± uj{kc) - Qc 



< Ct 



\kc + w\ 
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uniformly in w by 



sup 

P,1D,6 



L{k)dk 



< Ct . 



(6.22) 



6 — ^±{p,k) + irj \k + w\ 
This inequality is obtained exactly as (|6.21|) . 

Then we do the k^+i, ■ ■ ■ ,kn integrations using ( |6.1CI| ), eliminating all |-RjP factors in ( |6.20| ). 
This also collects 

t=^/2(ct)"-i(Clog*t)=^(Ct)"-"-2Ct(Ct)"-"t2H _ (ct)^-i/2(iog*i)3 _ 



Finally, we finish the estimate of ( 3.20 ) with the a,a,Pn integrations as in ( 6.14| ). This 
completes the proof of Lemma |6.2|. □ 



7 Fully expanded terms with big n and no recollision 



Here we prove ( p.30|) . The basic idea is that we consider a few disjoint crossing pairs (see 



Definition |6.1|) , which are well ordered (to make successive integration possible) and we plan 
to gain from each. 

The key is the right definition of "well-ordered" crossing pairs. The main difficulty lies 
in the fact that it is not true that "more crossings" give better bound. In particular, the 
contribution of the fully crossing pairing, 7r{j) = (n + 1) — j, is of order independently 
on the number of crossing pairs. 

So we need a more refined definition, which actually excludes the full crossing, and lets 
crossing pairs alternate with noncrossing ones. 

Definition 7.1 Fix a pairing tt G n„. A number a is calledpeak ifn^a—l) < 7r(a) > 7r(a+l) 
(Fig. 7) and it is called valley if it {a — 1) > 7r(a) < 7r(a + 1). The endpoints a = 1 and 
a = n are not considered peak or a valley. The valleys and peaks clearly alternate. We define 
a complete ordering on the set of peaks according to the value of their height 7r(a). 
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a-1 a a+1 



7r(a) 
Fig 7. A peak 

The number of peaks measures the complexity of the pairing. The proof of the following 
lemma is given in the Appendix p. 

Lemma 7.2 The number of pairings in n„ with no more than K peaks is at most ri^^^^{2K + 
2)". 

We want to gain a factor from each peak. The idea is that the factor i?7r(a)-i (together 
with Ra-i, Ra) will be eliminated by integrating out Pa~i,Pa at the expense of only log*t 
factors instead of the trivial estimate |-R7r(a)-i| < t. In order to estimate this integral, we 
make sure that only these three denominators depend on pa-i,Pa since we cannot estimate 
integrals with many denominators. This requires a certain combinatorial structure of the 
peaks. 

There is a difference whether u is constant or not. The u = (const.) case is simpler to 
integrate out and we will be able to gain a factor t^^ from each element of a monotonia peak 
sequence. 

In the uj 7^ (const.) case the propagators Rj,Rj have more complicated momentum de- 
pendence (see Appendix [A.2| for more details), hence we will need a stronger combinatorial 



structure for the peaks to be able to perform the integrations. The appropriate structure 



{"monotonia staircase") is introduced in the Appendix [A.l . 

We explain the proof of the uo = (const.) case in details. The necessary modifications for 
the nonconstant case are sketched in the Appendix |A.1|, |A.2. 



7.1 Combinatorics for the case of the constant lu 

The following theorem is a well known Ramsey-type theorem in combinatorics: 
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Lemma 7.3 Any sequence of different numbers of ajd + 1 elements either contains an in- 
creasing subsequence of {a + 1) elements or it contains a decreasing subsequence of {(3 + 1) 
elements. □ 



We will prove the following estimate. 

Proposition 7.4 Let n < N , m,ni & Ai{n, N). Suppose that u is constant and that vr G n„ 
has either an increasing sequence of k peaks or a decreasing sequence of k peaks. Then 



From these statements ( p.30| ) follows for the uj constant case. Choose k = 6, then all 
pairings which have a monotonic sequence of at least 6 peaks can be included into the second 
term in ( p.30[) . Choosing a = (3 = b and K = 5 • 5 + 1 = 26, it is clear from Lemmas |7]2 



and [T^ that with the exception of at most 54"7t,4-26+3 ^ (jn pairings, we are in the situation 
of Proposition |7^. The exceptional pairings are estimated by the apriori bound ( |5.1| ) to give 



the first term in ( |3.30| ). 



Remark: In general, we get t " with the exception of (Ck)^" pairings. 
7.2 Estimate for the constant uj case. 



Proof of Proposition Let ai < a2 < . . . < be the locations of the monotonic peak- 
sequence. We start with the expression ( |6.2|) . We estimate all the -Rj's trivially by Ct except 
with m = 1, 2, . . . , K, and except Rn- This gives a factor (Ct)""'^. In fact, with the 
help of we gain a factor from one of these estimates, at the expense of collecting 

{Pn)'^Y{j{kj)'^ factors. We also insert a factor 

1 <C" ((Pn)^ (7-1) 



this will help to secure a decay in (a). All these factors will be incorporated into the integration 
measures as before at the expense of changing M*{k) to L{k) = M*{k){k)^'^. 

Figs. 8-9 show the leftover i?'s represented by their electron momenta Pn{am)-i (bold lines) 
in case of an increasing and in a decreasing sequence of peaks, respectively. We want to gain 
a factor of t from each bold line. 
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p p p p p p 

Fig. 8. Increasing sequence of peaks 



i'a^-l 



Pa -2 Pa -1 



fa +1 



* a 

\ m 



3i(a J-1 7i(aJ 



;i(a^-l TKa^) 



P P 



Fig. 9. Decreasing sequence of peaks 



Now we express everything in terms of pj, j = 0,1, ... ,n. Recall that 

1 



R 



■7r(am)-l ~ 



a - e(p^(a^) +Pam-1 - Pam) + iV - ^7r(a^)-l 



where now 



3 ■ 7r(j)>T(am) 



i.e., it is independent of the electron momenta. Similar formula is valid for ^am- 

The important thing is that Pn{am) is independent of Paj-i and Pa^ for all j < m in the 
case of monotonically increasing peaks, and it is independent of Paj-i and pa^ for all j > m 
in the case of monotonically decreasing peaks. This can be seen from the expression 

j : 7r(j)>7r(am) 



and from the structure of the peaks. 
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Hence we can integrate out and pa^-i, m = 1,2, . . . , k, variables in the following order 



Va\ — l)Va\ ) Pa2 — liVa2 1 



if we have monotonically increasing peaks; and in the order 



PaK~li Pa^i Pai^-i — li Pa^-ii ■ ■ ■ ^ Pa2 — li Pa2i Pai — li Pai 



(7.2) 



if we have monotonically decreasing peaks. Notice that the dpa^ and dpa^-i integrations in- 
volve only the factors Ra,„-i, Ra^ and R-^^am)-! assuming that all other paj-i,Paj have already 
been integrated out for j < m in case of increasing peaks, and for j > m in case of decreasing 
peaks. 

The result of each integration is a C log* t factor by using the following inequality: 



L{pa^-2 - Pa^-l)L{pa^-l -PaJHpa^ - Pa^+l)dpa„,-ldpa„ 



{PaJ^ a - e{pa^_i) +17]- VLa^_i a - e{paj +ir]- VLa„ 

1 



(7.3) 



X 



< 



(Clog* tf{nuj) 



{(^){Pa^-2-Pa^+lY^^ ' 



The factor {pa.^)'"^ on the left hand side and the factor {a)~^ on the right hand side are present 
only for m = 1. The prototype of ( [7.3[ ) is the following estimate 

{p - u)-'^-^{u - vy^-^{v - q)-'^-^dudv {C\og*tf 



sup 

p,q,r,ei,§-^ {vy 01 - e{u) + if] 



6 — e(r + u — v) + iTj 



< 



2) 



(7.4) 



02 - e(f ) + if] 

with the choice p = Pa^-2, u = Pa^-i, v = Pa^, q = Pa^+i and r = p^(a,„). Again, if there is 
no factor on the left, then there is no factor (6*2)"^ on the right. Since 62 = a — fia™; 

we can estimate 

1 ^ C{niu) 



(02) - {a) 

We used that r = Pn{a,n) is independent of m = Pa^-i and v = pa„ so that we could take 
the supremum over r = p-n(am) outside of the integration. We also kept a momentum decay 
relation from the three L-factors using that L{k) < C {k)~'^~^ {k)~'^~^ and 



1 



< 



C 



— Pam-l) {Pa^-1 —Pa„,){Pa^ " Pa^+l) {Pa^-2 — Pa„,+l) 



Only part of the decay in L{k) was needed in 

For the proof of (|7.4]) we use (|6.19|) to perform the dv integration, then (|6.21|) for the du 
integration. 
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After having ehminated all remaining Rj factors, the rest can be done by successively 
integrating out 

PO,Pl, . . . ,Pa^_2,Par-l,Pai,Par+l ■ ■ • , Pa2-2, Paa-l , Paa , Paa+l • • • , Pn-1 

in this order (hat denotes missing variable). Notice that we saved a successive momentum 
decay for the remaining variables so we can use (see (|6.11|) ) 

sup / . ^ , . I < C log t 

c,q J |c — e[p) + trji 

with c = a — Qaj , P = Pj, 1 = Pj+i for all j ^ {am — '^,CLm — ^,cim '■ m = 1,2, . . . , k}, and 
simply q becomes for j = — 2. We collect (Clog* t)""^*^. 

At the end we are left with an integration identical to (|6.14 ). This completes the proof of 



Proposition 7.4. □ 



8 Amputated term without reabsorption 

Here we indicate the proof of ( |3.31| ) which is very similar to that of ( p.30| ). We explain the 
necessary modifications in the corresponding formulas leading to ( p. 301) . One has to check 
only Sections ^ and ^. 



Our starting point is Proposition Almost the same representation is valid for the 
amputated term TrD^^^Po V^ j^ except that in the set A4 {n, M) we additionally require 
mo = 0, /i(l) = 1, i.e., /o = and we also set Jq := 0. In particular the index is not part of 
the set J U J'^ and \J'^\ := n + \rn\ is reduced by one. This expresses the fact that there is no 
propagator 6**"'^^^°^ associated with the momentum pq. The definition of ~ ^ ( [4.21|) , ( [4.22|) 



is changed slightly: the products over j start from j = 1 and there is no Sq, Sq. 

The additional gain 1/t is due to the fact that the factors Rq and Rq are missing and that 
the set J*^ has smaller cardinality. 

In the proof of Lemma |5]^ one can easily see that there is one less factor of t"" in ( ^.5|) 
due to the decreased size of \J^\ in ( ^.4|) . The other t^~"- gain comes from the fact that the 
last integrand in ( |5.15| ) does not have a denominator, so the dki integration is bounded by a 
constant. 

In the proof of Propositions |7.4| the gain comes from the fact that Rq is missing hence it 
does not have to be estimated by Ct at the very beginning. Notice that Rq has indeed been 
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estimated by Ct since 7r(am) > 2 for all peaks. For nonconstant uj (see Appendix ^) the 
argument is similar: in the proof of Proposition |A.6| the first factor 5*0 is missing in ( |A.4| ). If 
77 (ai) = 1 and there are no Sj factors and (3 integration, then we gain from the missing Rq 
term. This completes the proof of (|3.31|) . □ 



9 One reabsorption 

In this section we prove ( |3.32| ). Recall that the measure 

p{n,N) -j-j definition of 

(see ( p.lO|) and (|2.13|) ) contains a double summation over measures /*(^'")ndfc,-, where 
m G A^a(r2, A^). Hence we can write 

TV 

KnU) = E E -^TNii) 

a=2 rn£Main,N) 



with 



*{m,a) , ^ , 



(9.1) 



After a Schwarz inequality we can symmetrize 



N 



ri,N 



a=2 meMa{n,N) \ 



V: 



m,a 

n,N 



using that the cardinality of A^a(n, A^) is bounded by . 
Therefore it is sufficient to show that 



lim sup Tie+ph ^^^Wro 



1 



2,\N 



(9.2) 



for any 2 < a < N and m G Aia{n, N). From now on we fix a and m and recall that 
= / \ {1, /i(a)}. Let i: := {1, . . . , AT} \ J,. 



We separate the measure of the reabsorpting phonon momenta in ( pTT 



dk 



#(m,a) 



beia 



n dkh) / dk,dk,^a){lldkhdkb+l)Fa{[k]j^,[k]l.: 



beJ 



with 



/ ( n dfc. 

be/a 



\{dh) n (}-x{h + h, 

b&Ia bjtb'ela 



(9.3) 
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and 



Fai[k]ia,[k]i-) ■= Xiki + k^(a))^ikTn) ]J [I - xikb + ki)j{l - xih + k^^^a))) , 

where in the notation we spht the variables k = {fci,...,/cAr} into two sets: {k} = [k]i^ U [k]ic 
and for any 5 C {1, 2, ... , A^} we let [k]s := {k^ : b E . We let 

r ^ 
J^{[k]i^,PN) := / dhdk^^a){U^^bdh+i)Fa{[k]i^,[k]ic)A{t,k,N;pN + Y.^J^PN) (9-4) 



be an operator in the phonon space, parametrized by [k]j^ and p^. Here A(t,k, N;pQ,pf^ 
is the kernel of the operator A{t,k,N) in Fourier space of Tie similarly to ( |4.29| ) and ( [4.30| ). 
Using the delta function 6(j)o — Pn — Yl!j=i kj^ we have 



AT 



N-1 



N 



-isj[e{pj)+Hph_ 



A{t,k,N;pr, + Y.k^,PN) :=A^/ ( n dpj)A*ip,k) f \ds,]^ {f[Q{k,)bk^e 



which is an operator acting on 7ip/i. Here we define A* by the relation A(p,fc) = A*(p, k)6(po — 



(9.5) 



PN-EUk,} (see (F27D)- 

The same formulas are valid for the other copy of T>^'^ in ( p.2| ); the variables are denoted 
by tilde. The result is 



#{m,a) 



{\{dk 



bGia 



#(rn.a) , _ ^ /• 

b] i ( n dkb) / dpN%{pN,PN] 



b€la 



X TiphJ" ( [A;] i,,Pn) iph ( [k] 



,P7V 



(9.6) 



We again used that the 5{po —po) coming from taking the electron trace can be replaced with 
^{Pn — Pn) and we integrated out p^. 

Similarly to Section |4.5| , we again notice that the phonon trace in ( ^.61 ) is zero unless there 
is a complete pairing between all the involved 2N momenta k, k. The pairing must respect the 
prepared immediate reabsorptions and the reabsorptions between ki, k^(^a) and ki, k^(^a), apart 
from an error term that is negligible in the thermodynamic limit (Lemma [4.4|) . Moreover, 
there is no pairing between kb and kb', b, b' G la (and the same for the variables with tilde) by 
the no-reabsorption condition built into the measure 

^#(m,a) Yi^j.^ (U). Hence the possible 
pairings are parametrized by a permutation vr G H^, where H^ is the set of all permutations on 
the set {2, . . . , a, . . . , n). The map tt* := on o : la ~^ la gives the pairing of the indices 
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of those A;'s and /c's that are not prescribed for immediate recoUisions similarly to Section ff3 
This means that a factor 



J2 Ilxl^Mi) -^M-(i))) 
can be freely inserted into ( |9.6|) modulo a negligible error. We obtain 

with 

( n ^^b) J [ n ^^b) J dpN%{pN,PN) 



(9.7) 



(9. 



b&Ia 



beia 



Peaks and valleys of vr G 11^ are defined exactly as before (Definition and Lemma |7]2 



remains valid if n is replaced with n — 2. The estimate (|9.2| ) follows from Lemma 7^ and the 
lemma below. 



Lemma 9.1 Let N > 3, a > 2, n > be integers and m G J^a{n,N) as before. For any 

(9.9) 



n gH'^ we have 



C^'Ni^;t)<l.iCXHf.h\og*tY 



Moreover, if n has k>2 peaks, then 



C^]^(vr;t) < f (CA^t)^ • -^(log*t)"+''+2 



(9.10) 



In case of ui = {const.), the power k — 2 can be improved to k. 



The proof of (|9.9|) and (|9.10|) are given in Sections |9.1| and respectively. These are the 



recoUision analogues of the apriori bound ( |5.1|) in Lemma ^]T] and Proposition |7J 



9.1 General bound for recoUision pairings 

For the proof of ( p.9|) , we start by symmetrizing the momenta in the definition 
Schwarz inequality within Tiph in ( |9.8|) 



TTphJ^[[k]j^,PN)lph y^[[k]i,,PN 



8[). Using a 



(9.11) 
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< 



Pn Iph 



J^[[k]la,PN)\ + ^T^'CphJ^[[k]la,PN)lph[j^[[k]la,PN 



We deal only with the first term, the second is identical. Hence 



#{m,a) 



beia 



#{rn,a) 



( n ^^b) / dpN%{pN,PN) 



bela 



~ \ 



3 = 2 



,Pn 



and we can freely integrate out all kf,, b G /„. Each k^ integration gives a factor so after 



these integrations 



1 



|A 



n-2 



#(rn,a) 



( n ^^b) / dpN%{pN,PN) TTphJ^([k]j^,pN)-fph[^{[k]l^,PN 



bela 



Notice that this bound is independent of vr. 

Now we write out JF explicitly using (p.4|) and (p.5|). We obtain 



X 



2N 



lA 



n-2 



#{rn,a) 



( n ^^b) / dpNle{pN,PN) / d/CidA;^(a)d/CidA;^(a) 



X 



X 



bela 

n dfcA+i) ( n dfc5dA;5+i)F,([A;],„, [^],c) 

6GJ feGJ 

AT-l „ AT-l AT 

n dp,)A*{p,k) / ( n dp,)A*(p,r) / Mf[d5,]f(nQ(^.)e 



-isje{pj) 



xTr, 



iV 



AT 



\{Q{k, 



Jsje{pj) 



J i=l 



with fc* := fcj if j G and /c* := kj if j E 1^, and we recall the definition of tj, fj from ( |4.34|) . 

Again, the pairing in the phonon trace must respect the recollisions prepared in the Fa 
factors and the prepared bk^bl. pairing for j G la (modulo an error 0(|A|^^)). These latter 
pairings yield a factor |A|"~^ from the n — 2 delta functions 6{kj — kj) = |A|. 

Now we proceed similarly to Section ^4.6| by integrating out the internal momenta kb, k^ 
for b ^ J, relabelling the external ones and expressing the relabelled electron momenta by 
the relabelled phonon momenta ( |4.24| ). We need a slight modification of the expression pj 
(j = 1, . . . , n — 1) compared to ( [4.24| ) 



Pj ■=Pn+ J2 ^rn , Pj ■= Pn + X(j < « - l)^a + ^ 
m=j+l 



(9.12) 



m— j + l 
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i.e., now these are considered functions of /c2, • • • , kn, ka and Pn- Similarly, we slightly modify 
the definition of fi^, j = 1, 2, . . . , n, (compare with ( |4.25| )): 

n n 

■= H (^m^{km) , ■= XU < a - l)^aUj{ka) + CTm^^l^m) • (9.13) 

m=j+l m=j+l 

The exact dependence on the a's are irrelevant. The important property is that Qj depends 

only on the momenta fc^+i, ■ ■ ■ ,kn and the same for Clj just ka is replaced by ka- 



The result of mimicking the argument in Section |4.6| leading to Proposition [4.3| is 

C^'Ni^;t)<a:^m (9.14) 

with 

C^Nit) : = A^^ J dpnle{Pn,Pn) J {f[M{k,,a,)dk,) (9.15) 



(T2,...,(T„;o-a j = 2 



X J M{ka,aa)dkaYrn,rn{t,Pn,Pn,k,k,aU{aa})+0{\A\ ^) 



where the definition of Y ( |4.22| )- (|4.23|) is slightly modified as follows. There is no permutation 



vr. The products over j in (|4.21| )- (|4.22|) start from j = I and there is no Sq, Sq. The sets of 
phonon momenta are k = [k]i^ U {ka} and k = [k]j^ U {ka}. The electron momenta p/s are 
functions of k as given in ( |9.12| ), similarly for the tilde variables. Finally, we use the definition 
of Qj, Qj (|9.13|) , in particular Y depends on {cr2, • • • , o"n, 5"a} which replace the variables a in 
the original definition of Y. 

The estimate of C^'^{t) is different for a > 3 and for a = 2; these cases will be discussed 
separately. The first case is similar to the crossing estimate in Section 3.3. [EY2]. The second 



case corresponds to the nested pairings. We recall the definition from Section 3.4 |EY2 



Definition 9.2 A pairing of the momenta {ki, . . . , /cat}, {ki, . . . , k]\f} is called nested if there 
exist indices ji < j2 < js < j4 such that {kj^, kj^) and {kj^, kj^) are paired, or the same is true 
for the tilde variables. Notice that this notion is defined in the original graph and not in its 
skeleton. 

In particular, a nested recoUision graph with a = 2 corresponds to ji := 1, J4 := yu(2) > 4 
and (j2;j3) with ^3 := j2 + 1 being the indices of any immediate recoUisions in between 
(1 < j2 < At (2) and j is even). 
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9.1.1 Pairing with a distant recollision: Case a > 3 



We assume that a > 3. In this case the recoUision pairing hne actuaUy crosses at least another 
pairing line. In | |EY2| | we treated this case together with the crossing estimates. The proof we 
give here is similar although in this paper we do not call it crossing pairing. 



To estimate Y in (|9.15|) , we bound all T's by a constant ([4.8| ) and use (|5.8| ). We obtain 
Cf'^{7i;t) < Cf^it) < (CA)2^t2|™|g^p f apMPn,Pn) r (9.16) 



— oo 
n 



X / I Y[M*{kj)dkj) M*{ka)dkaYl\R,\ \R 
i=2 j=i 



and we choose t] = t ^ in the definition of Rj , Rj ([4.23|) . 
For j > 1 we define 

n n 
Pj ■=Pn+ ^rn, ■= ^ra^^i^^r, 



m— j + l 



We start by integrating A;2, this involves Ri,Ri, explicitly (see ( |6.19|) ) 

M*{k2)dk2 



a - e{k2 + P2 + K) +ir]-Q2 - (^2^{k2) - Oa^^ika] 

1 



\ka ka\ 



a 



< 



Pa-1 



k, =-k 

1 a 



k =-k , 

1 a I 







k =k , 

a-1 a-1 


k =k ' k 

a+1 a+1 ' n 








' k 



N-l 




P P P^l P^ P 

12 ai a 

Fig. 10. One recollision after symmetrization and relabelling. 

Now we integrate out k^, . . . , ka-i in this order if a > 4. At the j-th step, 3 < j < a — 1, 
we eliminate Rj^i, Rj^i, using 

M*{kj)dkj 



a — e{kj + Pj + ka) + it] — flf — aju{kj) — OaUj{ka) 



(9.17) 
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and this gives (Ct)"~^ (this is vahd even if a = 3). This estimate follows from ( |6.10D after a 
Schwarz inequality. 

Now we do the dka, dka integrals to eliminate Ra-i and Ra-i'- 

M*{ka)dka 



I 



X 



1 



a - e{ka + Pa) +ir]- (JaUj{ka) - fit 

M*{ka)dka 



(9.18) 



{a — a) 



ka - ka\ a - e{ka + Pa) +ir] - aaUj{ka) - 



< 



The prototype of the dka integral is 

1 M*{k)dk 



sup 

c,q J \k — q\ 



< 



C\og*t 
{c-e{P)) 



(9.19) 



c-e{k + P) +ir]±uj{k) 

which is a straightforward strengthening of ( |6.21|) . The same estimate is valid without the 
\k — q\~^ factor, and that is the prototype of the dka integration. Finally, we use 



sup 



< 



C 



(9.20) 



P {c-e{P)){c-e{P)) - {c-c) 
to obtain the {a — a)~^ decay. 

Now we do dka+i, dka+2, ■ ■ ■ , dkn integrals; each gives a factor Ct and altogether we collect 
{Ct)"'~°'. In the j-th step (a + 1 < j < n) we eliminate Rj-i, Rj-i by using a Schwarz estimate 
and (lerrq ) as in ( |9T7| ). 

Finally we are left with 

1 r 7eiPn,Pn)dpn 



dada- v i , , , ~ 

oo {a — a) J la — e(p„) + zri\ \a — e[Pn) + tri\ 



<{C\og*t)' 



(9.21) 



Altogether we have 



(C log* 1)^01)"-^ {C log* t) (a)"-"(C log* t) 



(log* t )' 



icty 



for the integrals in (|9.16| ), which finishes the proof of (|9.9| ) in the a > 3 case. 



9.1.2 Pairing with a nested recollision: Case a = 2 



Here we prove ( |9.9| ) for the remaining case a = 2. In this case mi > 1 since we had a genuine 
recollision before relabelling. Therefore the original pairing was nested. 
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We start from ( |9.15D and perform explicitly the k2 integration. This involves the following 
A = Aia,p2,n2,a2) := / M(A;2, cxa) fi?i(a,pi, l^i, r/)!"'^' [T,(a - l]i,pi)l"'d/c2 (9.22) 



(9.23) 



using ([4.22|) and recalling that pi = P2 + ^2, ^1=^2 + o-2'-^(^2)- 
Lemma 9.3 With the notation above and for nii > 1, r] < 1 we have 

I A(a,P2, 1^2,^2)1 < , ^ , . 

(a - e{p2) - U2) 

Proof of Lemma \9.3i . We present the proof for uj = (const.) case here and explain the 
necessary modifications in Appendix |A.3 for the nonconstant case. The proof of the constant 
UJ case is a simple stationary phase calculation, while the nonconstant case is an integration 
by parts similar to Section 3.4 |[EY2 . 



For the constant u case we recall (|423| ), (gj) and we write 



J e-''<P'^ Tr,{a-n2-a2UJ,p 



M{p - p2, a2)dpds 



mi 



M{p - p2,a2)dp 



< 



after a change of variables. Using ( |4.8| ), ( |4.1(J| ) and that M has decaying derivatives up to 
order d we obtain by stationary phase argument (see (^4.13|) ) 

{sy/2{a-e{p2)-n2) 
After performing the ds integration, we obtain (|9.23| ) if c? > 3. □ 

The same estimate as (|9.23| ) is valid for the k2 integration which can be performed inde- 
pendently. We collected a factor (Ct)^™^^^ (with 77 := t~^) and we gain a (a — a)~^ factor, 
similarly to ( |9.2UD . The extra {p2)'^'^ factor can be compensated as before, using ( p.'S\) . 

Then we estimate all the other Rj and T factors in absolute value in ( p.l5| ), we use ( [4.8| ) 
and ( |5.8|) for j > 2, collecting (^(7t)2|nil-2m,i_ 'pj-^gn -^e follow the end of the proof in Section 
|9.1.1| by integrating dk^, . . . , dfc„ and further collecting (Ct)""^. Hence the total t power is 
N — 3, finishing the proof of ( |9.9D for the a = 2 case. □ 



9.2 Bound for recollision pairings with many peaks 



Here we prove (|9.10| ) of Lemma |9.1| . The argument is similar to the proof of Proposition 
in fact the recollision will not be used. 
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We start from ( |9.8D and follow the argument of Section ^?T| except the Schwarz inequality 
We obtain (compare with ( p.l5| )): 



j=2 



X / M{ka,aa)dka / M{ka,aa)dkaYjn,rn,n{t,Pn,Pn,k,k,aU {aa}) + 0{\A\ ^) , 



where the definition of Y (see ( |4.22| )) is modified as follows. The two products over j in ( [4.22| ) 
start from j = 1 and there is no sq, 5o- The electron momenta Pj,Pj, J > 1, are given in ( [4.24 ) 
and Qj, Qj are given in ([4.25|) with the modification that cr^-i(a) := cr^ in the definition of Qj. 
After the trivial estimates (^4.8|) and ( |5.8|) we obtain (compare with (|9.16| )) 



//•oo 
dpnle{Pn,Pn) / dada 
J —CO 

( n M*(A;,)5(% - Ku))dk,d~kj) / M*(A;,)dA;,M*(fc,)dL I] l^il l^il + Oi\A\-') 
where are in effect. 



This estimate is our starting point and it should be compared with ( |6.2D which was the 
starting point of the proof of Proposition [7.4| . 

For Lo = {const.) we follow the proof of Proposition |7.4| . Let ai < a2 < . . . < a,^ E la he 
the location of the monotonic peak-sequence as before. We estimate all |-Rj|'s trivially by Ct 
except -R7r(a™,)-i, m = 1,2, . . . , K. Of course the term Rq is missing, which gives the extra 
factor compared to the proof in Section ^72. 



The fact that ki + ka = ki + ka = does not change the dependence of pj and pj on 
k, k [j ^ Since Pq and pQ does not appear in Rj,Rj, j > 1, the only new momentum 
constraints due to recollision, po — pi + Pa-i — Pa = and Po — pi + Pa-i — = 0, play no 
role. Hence the argument in Section remains unchanged; the key point being that Pn{a,n) 
was independent of pa™-i,Pa™, in (|7.3|). Therefore the dpa^,dpa^-i integrations can be done 
successively; starting from ai or depending on the monotonicity of the peaks. 

For uj 7^ (const.) we follow the proof of Proposition A. 6 in Appendix A. 2 ; the proof is 
unchanged, just the dpo, dpo integrations and the corresponding .Ro, Ro factors are missing. 
□ 
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10 Computing the Wigner transform of the main term 



To identify the weak limit of W^main(^f-^{X,V), we test it against a function J E iS(R^ x 
R*^). We fix K > 5. We recall the definitions of -f~{t) (gJ^), Y^^^-, and 
<^e(^) = '^)- Then we have the following proposition which can be proven exactly 

as Proposition |4.3| . 



Proposition 10.1 For any fixed K > 5 and J E 5(R'^ x R'^) 



lim sup 

L— >oo 



T 



K-1 min{7V,A^} 
N,N=0 "=° 2ieA^(n,Af)^g_^(„^^) 7ren„ 



With N ■= {N + N)/2 and 



a 



,e{±} 

3 = 1, 



(10.1) 



(compare with (j^J^). The difference between C'j^- (^nd C^~^ from ( 4-1^ ) is twofold. We 
replaced the measure dh'n{Pn,Pn, k, k, S.) (see ^.2(\ )) with 



dUn{^,V, k,k, a) : = S{vo - V -J^^j) Je{^,VQ) ^e{v + ^, 



(10.2) 



It 

: Y[ ^i^j ~ kn{j))M{kj, aj)dkjdkj^ d^di 



and instead of l \4-24\ ) we considered pj , pj as functions ofv,C,,k,k as follows 



Pj = + 2+ J2 ' Pj = V - ^+ J2 ■ 

i=j+i e=j+i 



(10.3) 



With these notations, the definition of Rj,Rj ( \4-2^) and i\4-2^ ) are unchanged. □ 



The next lemma gives an estimate on all C* ~ (t), and a stronger bound if vr is crossing. 

° rn,rn,TT ^ ' ' ° ° 

Lemma 10.2 Let K > 5, n < N < K , n < N < K , m E M{n, N) , m E M{n, N), then for 
any tt G n„ 

(10.4) 



limsup|C*~ (t)|< ^^-^y 



where M := (A^ + n)/2, M := (A^ + n)/2 and < a < 1. Moreover, if n E Hn is a crossing 
pairing (n ^ id), then 



limsup|C;~ Jt)| < t~'/\C\Hf {log* ty . 

L—KX) 



(10.5) 
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Proof. For ( |10.4| ) we follow the proof of the analogous Lemma starting from ( |10.1D . 
We again split |y| = |y |i-'*/2_ ^j^e first factor is estimated in supremum norm using 
|r| < C^t^+*~7(M! M!) (see (U)). The estimate Q is modified as 

\C*^M < sup / dK{^,v,k,~k) 



Ml Ml. 



3=0 



j=0 



2-a 



with a modification of di/jr 

n 

n 

( n ^(^i - hu))M*{kj)dkjdkj) d^d 



(10.6) 



X 



The rest of the proof goes through until ( |5.2CI| ), just pn is replaced with f + f and p„ with 



|, i.e., (|10.3|) is used instead of ( |4.24|) , the t-powers are adjusted and d/i is redefined as 



dn{i,v,M) ■■= \Je{i,VQ) \ 6(vo-v 



\{l[M*{kj)dk,) d^di 



Let J*(0 •= sup^ \ Je{C,,w)\. After the d/c„ integration in the analogue of (p.20|) , we arrive at 



\ci^,M < (cx) 



N+N I 



M \ «/2 



Ml Ml 



t 



(|m| + |m|){l-a/2)+{l 



[ d^dv \%{v + iv-i) 



2+\N 



< 



d^dv - - ^) + %{v + + 1) 



(10.7) 



{Ml M!)"/2 

using \% {p,p')\ < lne{p,p)+%{p',p') \ by 7e > 0. Then ([10^) follows from (^JSj) and ( |L26D . 



The proof of (|10.5| ) requires very similar modifications along the proof of the analogous 
Lemma 16.21 The estimate (16. 21) is modified as 



C;~^(t)| < (CA)^+^tH+H sup / du:i^,v,kk) / da n 1^.1 / da J] 1^, 



j=0 



j=0 



keeping in mind (|10.3|) . 



It is easy to check that all estimates in the proof of Lemma |6.2| go through, since they 
were always valid uniformly in the incoming electron momenta denoted by p in ( |6.10 ), ( p. 11 ), 
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( |6.12| ), ( |6.21D and ( |6.22[ ), hence a shift ±| does not make a difference. Finally, instead of 
( |6.14| ) we have 



°° dada 
^oo (a) (a) 







{v+iy{v- 




a — e(v + 1^ 


+ ir] 




+ ir) 



< (clog*t)^ 

We again used ( |1.26|) and 

sup / j:(o(o'de = sup / sup \.u^,v)m'd^ < oo . 



e<l 



e<l J V 



This is an extension of ( |2.22| ), and it is clearly valid for J G iS(R'^ x R'^). □ 

From ( P.0.5D we know that only the direct pairing counts. To compute ~ we use the 
second formula for ~ (|4.22|) and that Vlj = Clj 

r _ /"OO 

= ^'"^ E UMi.vA.ka)e^'^' dae-*" (10.8) 



<Tje{±} 

j = l,...,n 



X 



j=0 



Introduce 



fi := V 



m=j+l 



j=0 



m = 1,2, . . . ,n — 1 



(10.9) 



and Vn '■= V. Note that pj = vj + |. We show that T^(q; — ^j,Pj) and T^{a — ^j,Pj) can be 
replaced with 



■^j ■.= To+{e{vj),Vj) , and j = To+{e{vj),Vj) 
modulo a negligible error (recall ( [4.11| )). Notice that \E'j depend on v, ^, fcj+i, . . . ,kn- 

Lemma 10.3 



with 



j = l,...,n 



(10.10) 



(10.11) 



0=0 



X / da e-**" n RT^'^^T' / da e**" H R?^'^' ■ 



j=0 
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We also have 



iimsup|c;*~,,(t)|< 

L— >oo 



2+\N 



, . (10.12) 

(M! M!)'*/2 

Proof of Lemma \L 0. Sj . The estimate (|10.12|) is proven exactly as (|10.4|) . For ( |10.10|) , we 
compute the difference C* — C**. We obtain a similar formula as ( |10.8|) , just Hj • O™^ is 
replaced with Y\j T {■ ■ ■)"^^ — Ylj "^J^^ and similarly for the factors with conjugate. We then 
estimate this expression exactly as in the proof of (|10.4|) with a = 0. 

Let 6^ : R ^ R be defined as 9{s) := \s\ for |s| < 1 and 9{s) = 1 for |s| > 1. Using ( |1.15|) , 



]8D and (|4.9| ), we see by a telescopic estimate that 



j=0 



j=0 



j=0 



9(a-n,-eip,))+Om{P3) + {0) 



(10.13) 



C 



For the first term in (|10.13|) we notice that 

in other words, the net effect of a factor 6(^a — Qj — e{pj)^ is that it neutralizes the singularity 
of an \Rj\ factor and still keeping its decay property. Hence, effectively, rrij is decreased by 
one (see (|5.8|) ). We gain a factor t from this decrease (notice the dependence of the t power 
on |m| in ( |10.7|) ), but we lose r]"^/^ = t^/^ in (|10.13|) . This gives the extra factor t"^/^ in the 
error term of (|10.10|) relative to the robust estimate (|10.4|) . The special case rrij = require 
small modifications which we leave to the reader. 

For the second term in (|10.13|) we use that / |^|(0'^e(Od^ = 0(e) to gain r^'^/^e = t"^/^. 
The factor (pj) can be absorbed into the decay of M* and 7e as before. □. 



Combining Proposition |10.1|, Lemma |10.2| and Lemma |10.3| we obtain 



Proposition 10.4 For K > 5, 



lim sup lim sup 

e^O I/— >oo 



T 



K~l min{N,N} 

E E 



N,N=0 



n=0 



E E 

rneM(n,N} rn€M{n,N) 



. 



□ 



Using the factorials in the estimate ( |10.12| ), we can extend the summations over all m, m; 
the error term goes to zero as K ^ oo: 



lim sup lim sup lim sup 

K—>oo £— >0 L— >oo 



T 



n=0 



0, 
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where 



C{n,t):= E E C^^Jt). 

mo,...,m„=0 mo,...,mn=0 



The summations over mj,mj give a geometric series in (|10.11|) . We obtain a shift in the 
denominators of Rj {"one loop renormalization"): 

1 



a — e{pj) — Vtj — X^'^ j + ii] 



;io.i4) 



We can change back the da and da integrations into subsequent time integrations since \l/'s 
are independent of them (see the similar identity in ( [4.50| )). We obtain 



„^^s^x J -Jo Jo 



;io.i5) 



j — l,...,n 

X exp \ — i 



The rest is similar to the end of Section 4 in ||EY2|| and we will skip a few technical details. 



We introduce the new time variables aj := {sj + Sj)/2, bj := {sj — Sj)/2 for j = 0, 1, . . . , n 
and we have 

n— 1 



2) 2 



/It p I ii ±. 

( n dv,)di Ui. ^O) / ( n ^(^.- - ^^+1' ^3+l))le{^n + 

j = l,...,n 

X / ( n - E exp - ■ 5: a, Ve(t;,) + 2 a, (A^M,- + 0(e')) 

° 1=0 .7=0 I .7=0 .7=0 J 



n n 



i=0 j=0 j=0 I ^ 

Here we used pj = "^j + |, = f j — f ^iiid the second order Taylor expansion of e{y) 



with a uniform error bound (see ( 1.15 )). 

We rescale all microscopic variables into macroscopic ones, i.e., aj :- 
t = e~^T and recall that e = A^. We define 

/ n-l \ n-l 

Xa{b) ■.= x\ - oo < E ^ «o n Xi-aj < bj < aj) 



3=0 



3=0 
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for any a := (oq, ai, . . . , a„) and b := {bo, bi, . . . , bn-i)- Then 

n-l 



C{n,t) 



,e{±} i=0 i=0 

/•T " " /-OO 



n-l 



j=0 [ 

n 



X 



exp - < ■ E«iVe(t;,) + 2^«,-(m,- + 0{eC)) W^M,Vn) • 

I j=0 j=0 ) 

It is easy to see that the error terms are neghgible as e — > since T is fixed and C, is essentially 
bounded by the decay of J. The dbj integrations give 2Tx6{e{vj) — e(f„) + Vtj^ delta functions 
on the energy shell as e — * 0. 

Following the analogous calculations on pp. 709-710 of [[1l)Y2|| , using the macroscopic profile 
of the initial state ( |1.25| ) and recalling the definition of the Boltzmann collision kernel ( |1.27| ) 



a{y, U) := 27r ^ MiV - U, a)6(e{V) - e{U) + auj{V - U 
we obtain {Vj = Vj) 



a=± 



lim lim C(n, e~^T) = [ dXdVo . . . dK J(X, Vq) ( ft 

e^O L^oo J Jo \ 



daA5iT-Y: 



j=0 



■ n-l 



10.16) 



■ j=o 



j=0 



We also rearranged the energy conservation delta functions as 

n— 1 n— 1 

j=o j=o 



Noticing that 



-21m^j = J a{U, Vj)dU 



by ( |4.12| ), we see that ( |10.16| ) is exactly the n-th order term in the Dyson series solution of 
(|1.2|). This completes the proof of the Main Theorem. □ 
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A Case of general u 



Here we show how to modify the proof of Lemma if uj is not constant. The assumption 
that UJ = (const.) was used only in Sections ^ and |^. 

The proof of ( p.30|) for u ^ (const.) will be given first. The necessary combinatorial 
structure is introduced in Section |A.1| and the key analytic estimate (Proposition |A.6| ) is 
proven in Section [A.2| . These sections show how to modify the argument in Section |^. The 



modifications for the proof of ( |3.31| ) were already explained in Section |^ both for constant and 
nonconstant u. 

Most of the proof of (p.32|) in Section ^ is valid for arbitrary u. The assumption u = 
(const.) was used only when we estimated the nested pairing in Section |9.1.2| (proof of Lemma 
[9.31 ). The technical modifications for nonconstant u are given in Appendix [A.3| . 



A.l Combinatorics for general cu 

Since the momentum-dependence structure is more complicated, in order to perform the 
succesive integration of Section |^ we need more control on the structure of the pairing. 

Definition A.l Fix a pairing vr G Hn- A sequence of consecutive numbers a,a + l,a + 
2, . . . ,a + h G {1,2, ... ,n} is called a down-stair of length h if 7r(a) > 7r(a +!)>...> 
7r(a + h); and it is called an up-stair of length h if 7T(a) < 7r(a + 1) < ... < 7T(a + h). We 
always assume that h > 1. 

Notice that the elements in a stair must be consecutive, it is not sufficient if they just form a 
monotonic subsequence. 

Definition A. 2 A set of k up-stairs; 

ai, ai + 1, . . . , ai + hi ; 

02, a2 + 1, . . . , a2 + ^2 ; 



is called an increasing ^-staircase if it satisfies the following properties: 

66 



(i) ttj + hj < ttj+i for j = 1,2, ... ,K - 1; 

(ii) Each aj + hj is a peak; 

(in) The t: -images of the peaks are increasing: 

7r(ai + hi) < 7r(a2 + /i2) < • • • < 7r(a« + Hk) ; 

(iv) There is no peak p with Uj < p < aj+i, 7r{p) > 7r(aj + hj), j — 1,2, . . . , k — 1; 

(v) The images of the stairs minimally overlap, i.e., 

7r(aj+i + 1) > 7r(aj- + hj) > 7r(aj+i) j ^ 1,2, . . . , k — 1 . 

The numbers ai + hi, a2 + h2, ■ ■ ■ , + h^ are called the tips of the stairs, the numbers 
oi, 02, . . . , O/t are the bottoms of the stairs. 
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Fig 11. Graph of a permutation with an increasing K=3 staircase. 



In other words, an increasing K-staircase consists of k, up-stairs with lengths /ii, /i2, ■ ■ ■ , 
in such a way that every up-stair ends just above where the next one starts, i.e., if we removed 
the bottom from each stair then the stairs would not overlap. (See Fig. 11). Moreover, there 
is no peak between two stairs which would be higher than the tip of the lower stair. This 
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implies, in particular, that there is no dot at all in the shaded regime, i.e., there is no p with 
ttj < p < Oj+i, 7r(p) > T:{aj + hj), j = 1,2, . . . , k — 1. Note that the bottoms are not necessarily 
valleys. 

Similarly we define: 

Definition A. 3 A set of k down-stairs; 

ai,ai + 1, . . . ,ai + hi ; 

a2, 02 + 1, . . . , a2 + /i2 ; 



Ok, Ok + 1, . . . , Ok + K 

is called a decreasing K-staircase if it satisfies the following properties: 
(i) Uj + hj < Uj+i for j ^1,2,..., K- 1; 
(a) Each ttj is a peak; 

(Hi) The t: -images of the peaks are decreasing: 

7r(ai) > 77(02) > . . . > 7r(a„) ; 

(iv) There is no peak p with aj < p < 7r(p) > 7r(aj+i), j = 1,2, . . . , k — 1; 

(v) The images of the stairs minimally overlap, i.e., 

n{aj + hj — 1) > 7r(aj+i) > 7r(aj + hj) , j = 1,2, . . . , k — 1 . 
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Fig 12. Graph of a permutation with a decreasing K=3 staircase. 



The following Proposition shows that essentially every permutation has either an increasing 
or a decreasing K-staircase. This is again a Ramsey type theorem. However, the estimate is 
very bad in k, so practically it can be used for finite k only. 

Proposition A. 4 For any fixed k, with the exception of 

^4.4-1+3^2 -4"-^ + 2)" 

pairings, a pairing it G Hn either contains an increasing or a decreasing n-staircase (or both). 

Proof. We need the following lemma whose proof is given in the Appendix 

Lemma A. 5 Suppose that a permutation has at least ("^'^j peaks. Then it has either an 
increasing {a + l)-staircase or it has a decreasing {(3 + l)-staircase (or both). 



From this Lemma and Lemma |7]^ the proof of Proposition |A.4| follows easily. Choose 
a = /3 = K - 1, and use that (^^Zi) < 4^^"^ □ 



In the next section we will prove 
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Proposition A. 6 Letn < N, m,m G A4{n,N). Letuj ^ {const.), satisfying the assumptions 



in Section \1.2^ and assume that vr G n„ has either an increasing K-staircase or a decreasing 
K- staircase. Then 

\c^,^At)\ <t"'^" (cxHf {log* tr^^^' . 

From these statements (|3.30|) follows for the u ^ (const.) case. Choose k, = 8, then all 
contributions which have a monotonic 8-staircase can be included into the second term in 
(g). The exceptional n^'^ +^{2 ■ A'^ + 2)'' < pairings can be estimated by the apriori 
bound from Lemma |5.1| and be included in the first term in (|3.30| ) . 

Remark: In general, we gain t"*^"*"^ with the exception of ra'^'^'' ^^^(2 ■ A'^^'^ + 2)" pairings. 



A. 2 Estimate of the indirect term for the nonconstant uj case. 



Proof of Proposition \A.(^ . The difficulty is that there are momentum dependences within 
the argument of u. This makes it impossible to eliminate all the designated R factors by 
integrating out the peak variables as in ([7.2|). In Section |7]^ when we integrated out Pa„-i 
and pa^ , we used that only three denominators (two R and one R factors) depended on these 
variables (see (|7.3| )). In case of nonconstant u, still only these three denominators depend on 
Pam-i Pam the electron kinetic energy is concerned, but due to the cumulation of 

Lj's in the form of ±^7(^^+1 — ± — Pi+s) ... in Re, all Rj factors depend on Pa^-i 

and pa^ for j < am- 

Hence we have to integrate out all pj in consecutive order: pi,P2, ■ ■ ■■ This means that the 
designated -R7r(a,„)-i factors are eliminated in an alternating order with the regular Rj factors. 
One has to ensure that when we integrate out pj with some flm-i < j < — 1, the designated 
i > m, denominators are not affected since we cannot integrate out a variable which 
appear in many R factors. 



The staircase construction of Section |A.1| is designed to ensure these restricted dependences 
so that the pi,p2, . . . integrations could be done in this order. 

There is one additional difficulty: the staircase controls the ordering of the momenta 
between the first and the last stairs, but it does not control the ordering before and after the 
staircase. The natural idea to consider the "first" staircase does not work (the first staircase 
could be too short). 
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So we have to separate the staircase from the rest, unless ai = 1 or 7r(ai) = 1. We use the 

in. 

jJo- 



semigroup property of the measure /**[ds ^"- 



t* rt 



in+m 



(A.i: 



and we will lose a t-factor by estimating the outer integral trivially. This is why the estimate 



in Proposition |A.6| is weaker by a factor than that in Proposition 



A. 2.1 Estimate for increasing staircase 

Assume that vr has an increasing K-staircase with aj bottoms and aj + hj tips for 1 < j < n 
(see Definition |A.2D . Our starting point is ( [4.44| ). We use (|A.1|) for the time integrals in ( [4. 44)) 
as follows: 



t* 



N 



j=ObeijUi^ 



t I rs* 

ds[ / [ds^; 
\ Jo 



fi{ai)-l 




ai — 1 



Q-'i-Sb[e{p,i(j)+Xbkb)+^j-^ai-i+Xb^b^^{kb)] 

L j=0 6G/,U/| 



X 



(t-s)* 



J~J J~J Q-iSb[e{Pf,(j)+Xbkb)+^]+Xbf^b'^{kb)] 
lj=aib€ljUl^ 



(A.2) 



i.e., we separated the dsj, j < /i(ai), integrations from the rest by prescribing their total 
sum to be s. We also pulled out a factor e~*'^^°i^^ which shifted all fij's in the first group of 
propagators by Qai — 1- Notice that 



ai — 1 



m=j+l 



(A.3) 



depends only on A;j+i, . . . , ka^_i, and Vtj for j > ai depends only on ka-^+l, . . . ,kn. 

The decomposition of the dsj integrals is similar; we separate the dsj, j < /i(7r(ai)), inte- 
grations (their sum is s) and pull out a factor e**^'^(''i)-^ We again see that the /c- dependences 
of 



f2j — ^TT{ai)-l 



E 

m 

j<7v{7n)<7v{ai) 



for j < 71 [ai) — 1, and those of 



= ajnuj{km) 

m : TT{m,)>j 
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for j > 7r(ai), are separated. 

Now we use ( [4.50| ) for each time integral in ( |A.2| ) separately, we integrate out kh, b E J, 



and similarly for ki,h E J as in ( [4.51| ). We then estimate everything by absolute value, we 



use (|4.8| ) and (|5.8| ) and we estimate the two outside time integrals (ds, ds) trivially. In this 
last step we lose an extra . We also integrate out p„ and all kj and express everything as a 
function of kj and Pn- The result, similarly to ( |6.2D , is 



\C^,^nM < {CXr''t'+'\^\snp [diy:{p^,p^,k,k) (A.4) 

/oo ai-1 n _7r(ai)-l _ " _ 

dad/? n I'S'jl n l^il / dad/3 J] I'^'il H l^il ' 



°° j=0 j=ai - ^ i=0 i=7r(ai) 

where in addition to recalling the definition of Rj ( [4.23| ), we define 
S'j := S'j(/3,pj,nj - ^^1-1,7]) : = 

Sj := Sj{(3,pj,nj - l]^(a^)_i,?7) : = 



1 



(3 - e{pj) - {ilj - fi^(ai)-i) - 
and, as usual, Pj,Pj are viewed as functions of Pn and k ( ^.24| ). 



If Oi = 1 or 7r(ai) = 1 then there is no need for separation, i.e., there are no additional d/3 
or d/3 integrations. In this case (|6l^ ) can be used directly. We will not discuss this simpler 
case in detail. 

We introduce 

bj := aj + hj , 1 < J < fi^ 
for the location of the tips. All factors in ( |A.4| ) with tilde except Rn and j = 



1, 2, . . . , K, are estimated trivially by (Ct), this gives (Ct)""'^. Moreover, we make sure that 
we gain a {a)~^ {(3)'"^ factor from these L°°-estimates, exactly as in the constant u case, using 
( |6.16| ). Notice that we need to gain a second "^-factor to make the d/3 integration finite. 



This is possible only if 7r(ai) > 1, but there is no need for djS integration at all if 7r(ai) = 1. 
We also insert an explicit extra Ylj{kj)~^ decay. Since all (3 denominators are estimated in 
norm, at the end we can perform the d(3 integration. We get 
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with 

L*{k) ■= M*{k){ky'^+^^ (A.5) 

so L*{k) < C{k)^'^'^^^. Notice that we also inserted (Pai-2)~^(Pai-i)"''"^(P6i-i)~^ using ( |6l8D 
to prepare for the decays in /? and a at the expense of increasing the (j9„)-power and using up 
a few Y\j{kj) power from M*. If oi = 1, then there is no need for insertion. 
We express the electron momentum in Rn{b„,)-i as Pn{brr,)-i = Pnibm) + Pb,n-i - Pbm- 




P P P P 

Fig. 13. Increasing staircase. 



Fig. 13 shows an increasing staircase; the bold lines indicate the electron momenta in 
those Rj's which we kept (j = 7r(6,,„) — 1, m = 1, . . . , k). 

Now we change variables; we use po,Pi, ■ ■ ■ ,Pn-i instead of ki, k2, ■ ■ ■ , kn- Notice that for 
each ni < K the factor _R^(fc^)_i, i.e., the momentum p-K{bm)-i and do not depend on 

Pj, di ^ j < bm — ^- This follows from the staircase construction. 

First we estimate the factor L*{pa^^i — Pai) by a constant. In this way the p-dependence 
of the remaining L* factors are separated. The inserted (Pai-i) '^^^ factor ensures that the 
integration of the momenta before the staircase is still finite. 

We then integrate out pa^,paj^j^i, . . . ,Pbi-2, these variables do not appear in Rn{bm)-i, = 
1, . . . , K, or in Sj, j < cti — 1 (see ( |A.3|) ). We get a (Clog* tY^~^ factor by ( |6.11| ) and we have 
eliminated Rj, j = ai, ... ,61 — 2. For the Pbi-i and pb^ integrations we are in a situation 
similar to (|7.3 )- We explicitly write out the factors Rb^-i, Rbi and -R,r(iii)-i that are involved 
and we obtain 

L*{Pbi-i -pb^)L*{pb^ -Pbi+i)dpfei_idp6, 



{pbj-iy a - e(p6,_i) + ir]± uj{pb^-i - PbJ ± (^{Pb^ - Pb^+i) - ^b^+i 



(A.6) 
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X 



1 



a - e(p^(fei) + Pb^-i - Pbi) +ir]± uj{pbi-i - Pbi) - ^7T{bi)-i ~ ^^i+i) 



X 



< 



Notice that p-K{b^) is independent of pa^.pa^^i, ■ ■ ■ ,Pbi, i-e., of all the variables we have inte- 
grated out so far. The prototype of this inequality is 



sup 

g,r,e {vy 9 — e{v) + 17] ±. uj{v — m) ± uj{u — q) 

1 1 



(A.7) 



{C\og* tf 

(d) 



9 — e{u) + i?7 ± uj{u — q) 9 — e(r + v — u) + irj ± uoiv — u) 



< 



(with V = Pb^-i, u = pb^, q = Pbi+i, r = Pn(bi), 9 = a-Qbi+i, 9 = a-Q^(^bi)-i) which is proven 
similarly to (|7.4| ). First we perform the df integration using (|6.19 ), then ( 6.21|) estimates the 
du integration. We can change the decaying factor {a — Qbi+i)^^ in (|A.6|) into (a)^^ at the 
expense of the product of (kj)"^ = {pj — Pj-i)"^ factors (j = 6i — 1, . . . , n) similarly to (|6.13|) . 

Next we integrate out pbj^+i,Pbi+2, ■ ■ ■ ,Pa2y ■ ■ ■ ,Pa2+h2-2 = Pb2-2- By the staircase property 
(especially (iv) in Definition [A.2| ), these variables do not appear in the remaining factors 
Rnibj)-!, j > 2, hence they freely give C71og*t factors and they eliminate Rj, j = bi + 
1, . . . , 62 — 2. Now we integrate out ^62-1 Pfe exactly as in (|A.6| ), etc. 

Once we are done with all Pb^-i Pbm < rn < k) integrations, then there are 
no more R factors left. So we integrate out the remaining variables, po,Pi, • • • ^Pai-i, and 
Pb^+i, ■ ■ ■ ,Pn-i in this order; it is easy to see that at each step only one Sj or Rj factors 
depends on these variables. At the Pai-2 and Pai-i integration we use 



L*{pa^^2 -Pai-l)dpai- 



(3 - e(pai-2) +ir]± Uj{pa^-2 - Pai-l) (Pai-2) 



< 



Clog*t 



and 



/ 



< 



Clog*t 



(A.8) 



(A.9) 



P-eiPa,-l)+tV {Pa,-lY+^ 

that follow from ( |6.12| ) and ( |6.11|) . The extra factor makes the d(3 integration convergent. 

Finally from the a, a and p„ integration we obtain a factor (Clog*t)^ similarly to ( |6.14|) . 
This completes the proof of Proposition [A.6| for the case of increasing K-staircase. 
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A. 2. 2 Estimate for decreasing staircase 



The proof for the decreasing staircase is very similar. Here we chop off the end of the expansion. 
Let bj := aj + hj as before, but now these are the bottoms of the stairs. Then instead of ( |A.2[ ) 
we split the Sj variables into two groups: j = 0, 1, . . . , — 1 and j = fi{b,^), . . . ,N. The 
result, analogous to ( |A.4| ), is 



dad/3 n l-^jl n l^i 



j=0 



dad/3 n i-^ii n 1^. 



■3\ ' 



j=0 



j=Tr(b^) 



where Sj and Sj have been redefined accordingly: 

1 



Sr--- 



Notice that we pulled out the factor e"**^° and e**^'^('"')-i in the analogue of ( |A.2D . In this 



way, 5*^ depends on po, ■ ■ ■ ,Pj and 5*^ depends on pj, . . . 

All factors with tilde except i?„ and >S'7r(a™,)-ij m = 1,2, . . . , k, are estimated trivially, and 
we can ensure the decaying (a)^^(/3)^^ terms by inserting approriate factors at the expense 
of {pn) and {kj) powers as before. We get 

' -Q L*{kj) S{kj - ~k^(^,))dkjdk,^ 



{k,Y 



dpn{PnY'^^^'^le{Pn,Pn) 



n n i^ii n \s^{am)-\ 



in=l 



{Pa,-2r+HPb.-iy+'{Pa.-l)'{PbJ^ V ,=0 / V,=.. 

and we will again use the variables Po,Pi, ■ ■ ■ ,Pn-i instead of ki, . . . , k^. 

We first estimate the factors L*(pa^_2 — Pai-i) and L*(pb^„i —pbj by a constant. Then the 
integrations of Pb^-i, ■ ■ ■ ,Pa^+i can be done successively in this order; no 5*^ factor depends on 
these and we collect a {C log* t) factor from eliminating each Sj. We then perform dpa^-i, dpa^ 
integrations, very similarly to ( [A. 61 ): 

L*{Pa^-2 - Pa^-l)L*{Pa^-l ~PaJdpa^-ldpa, 



{Pa^-lV P - e(Pa«-l) +ir]± Uj{Pa^+l -PaJ± Uj{Pa^ - Pa^-l) " fia«+] 

1 



(AlO) 



P - eiPaJ +iV± ^(Pa.+l -PaJ- ^a, 
1 



X 



a 



- e{pn{a^) + Pa^-l " Pa«) + If] ± Uj{pa^ - Pa, 



Tr(b^)-l 



< 
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Again, we used that pn{a^) does not depend on the momenta integrated out so far. We perform 
Pa„_2,Pa^-3, • • •, Pb^-i, ■ ■ ■ ,Pa^-i+i successively; no remaining Sj factors depend on them. We 
then perform Pa^_i,Pa«_i-i exactly as in ( |A.10| ), etc. Once we are done with pai,pai-i, there 
are no Sj factors left, and the rest can be integrated out successively: po,pi, . . . ,Pn-i- The 
dpai-2 integration uses the (dpa^_2) '^"^ decay since there is no L* factor left. We also gain a 
{l3)~^ factor from this integration similarly to ( |A.9| ). The decay {jS — Qa^+i)^^ is changed to 
using the extra decaying factors similarly to ( |6.13| ). This ensures the finiteness of the /3 



integration. From the pb^ integration we gain an extra (a) ^ because of the inserted (pb^) ^• 
Finally the a,a,Pn integrations are done exactly as in ( |6.14| ), which completes the proof 



of Proposition |A.6| . □ 



A. 3 Estimate of the nested term for the nonconstant u case. 

In Section |9.1| we proved ( |9.9| ) for a > 3 and for a = 2 with u = (const.). Here we outline 



the argument for the remaining case when oo ^ (const.), a = 2. Only the proof of Lemma |9l3 



needs to be modified, the rest of the argument from Section ^TT] is unchanged. 

Proof of Lemma \9. for uj ^ (const.). The proof is an extension of the argument in Section 
3.4 of ||EY2|| and we only indicate the main steps. We write 



A = i"'+^ J M(k) s'^e''^''-'^'^P'''^+''''>[Tr,(a - u(k),p + k)]'^dsdk , (A.ll) 

where we omitted a2 and we let m = mi, p = P2, and ^2 = for simplicity. By a partition 
of unity we divide the space into cubes Q of size g. Correspondingly, we can replace M(k) in 
( |A.11| ) by MQ(k) that is supported on a cube Q and M = J2q Mq. We can assume that Mq 
is as smooth as M. Moreover, ||Mq||oo < C(dist(0, Q))"^^ and similar estimates are valid for 
its derivatives, since Mq inherits the size of M on Q. We denote the corresponding expression 
by 



Aq = AQ(a,p) := J MQ(k) s-e^K"-*(M)+^'7) [t^(c, _ uj(k),p + k) 



m 



dsdA; 



and we neglect the p variable which we consider fixed. 

After a diffeomorphic change of the k variable we arrive at one of the following normal 
forms depending on whether Q contains a critical point or not: 

Aq = j EQ(a, k) s^e^^^^-^'+'^MsdA; ; (A.12) 
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or 



Aq = J EQ{a,k) s'^e^^^^-^-^+^^MsdA; , 
where m = is a constant unit vector, a = ag G R, and we defined 

EQ{a, k) := I'^+'MQiijik)) [T,(a - ^(^(fc)), ^(A;) + p)] "j{k) . 

Here the function expresses the change of variable, J is the corresponding Jacobian. These 
are as smooth functions as $ and uo and they depend on Q and p with uniform bounds. 

We will discuss the more complicated first case, the other case is similar. We rewrite the 
dfc integration as follows 

j e-'''''EQ{&,k)dk = (27r)-'^/2^°°g-is^|^^^ ^EQ{a,k)da{k)^^'^/^-^d^ , 
where da{k) is the normalized surface measure on the sphere {k : k"^ = i)}. Let 

HQia,k): = J^^^^^EQ{a,^daCk)] (A.13) 



[l-l)k 



-2 



d 



EQ{a,k)da{k) - i— 



EQ{a, k)da{k) 



Using integration by parts 



e-''^^EQ{a,k)dk = s-^ I e~'"'^HQ{a,k)dk 



hence after undoing the ds integration in ( |A.12[ ) 



HqIo, k)dk 



[a — k"^ + irj)^ 
We show that 

\HQia,k)\<CQir^-'/' + \kn. (A.14) 

Hq consists of two terms ( |A.13|) , the first one is easy since EQ{a,k) is uniformly bounded 
from ( |4.8|) . For the second term, we use the bound on the derivative of T^(q;,p) in p and a 
( [4.91) . All the other factors in Eq have derivatives bounded by a constant. 

Since k i— >• HQ(a,k) is compactly supported, ( |A.14| ) gives an estimate of order Cgt™"^/^ 
for Aq (with 7] := t~^), where the constant Cq depends on the cube Q and it behaves as 
Cq < C(dist(0, Q))"^"^. In particular, these estimates are summable over all the cubes Q of 
the partition of unity. This will complete the proof of Lemma p.3| . □ 
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B Proofs of the combinatorial Lemmas 



Proof of Lemma \7.^. We consider a pairing vr G n„ that has exactly K peaks at the locations 



oi, a2, . . . , ttK- Such pairing has K — 1, K or K + 1 valleys at 60, ■ ■ ■ ,bK such that 

bo < ai < bi < a2 < ■ . ■ < bx-i < clk < bx , 

where 60 and bx might not be present. 

The number of possible peak-and-valley configurations, including their heights is at most 
^4K+2^ since one only has to prescribe the values aj,7i{aj) and bj,TT{bj). 

Once the peak-and-valleys are fixed, we consider the set 

5 = {l,2,...,n}\U(vr(a,)U7r(6,)) . 

j 

We define the following sequences: 

D,: = (7r(a, + l),7r(a,+2),...,7r(6, -1)) 
U,: = (7r(6, + l),7r(6, + 2),...,7r(a,+i-l)) 

for j = 0, 1, . . . , -f^ and ao = 0, ax+i = n + 1 for definiteness. Clearly S = Uj (Uj U Dj^ is a 
partition; and Uj (Dj) is a monotonically increasing (decreasing) subsequence. This partition 
consists of at most 2K + 2 non-empty sets, we label them with distinct labels. For any 
permutation tt with fixed peak-and-valleys, we assign to every element in S the label of its set 
in the partition. This can be done by at most (2K + 2)^^^ < (2K + 2)" ways. Notice that once 
the peak-and-valleys and this assigment are fixed, the permutation is uniquely determined, 
since the exact order within Uj and Dj are determined by the monotonicity. Hence the number 
of such permutations is not more than n^^^^{2K + 2)". If we consider pairings with at most 
K peaks, then we have to add these numbers up for K = 1,2, . . ., which gives a total number 
at most n^^+^{2K + 2)" □ 



Proof of Lemma \A.^ Let f{a,j3) be the smallest number such that any permutation with 



at least /(a, (3) peaks has either an increasing (a + l)-staircase or a decreasing (/3-|-l)-staircase. 
We show the recursion 

/(«,/? + 1) < /(« - + 1) + f{a,P) (B.l) 
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and the relations 

/(«,!) = « + /(I,/?) = /? + !. (B.2) 

From these recursive relations we easily obtain the estimate f{a,(3) < We first notice 

the following fact: 

Observation: If we have an increasing K-staircase with tips at oi < 02 < . . . < Ok, and if 
there is a peak a > which is higher than the highest tip, i.e., Tcla^) < 7r(a), then we also 
have an increasing (k + l)-staircase. Similarly, if we have a decreasing K-staircase with tips 
at ai > a2 > . . . > a^, and if there is a peak a < ai which is higher than highest tip, i.e., 
7r(ai) < 7r(a), then we also have an decreasing (/t + l)-staircase. 

For the proof of this observation in the increasing case, suppose that we have a peak a > 
with Ti{af^) < 7r(a). Then, we simply choose the smallest number h > a^, such that 6 is a peak 
and 7r(6) > n{a^). Let b be the tip of the (k + l)-th up-stair and let < c < b be the biggest 
number such that tt{c) < Tria^). Such number exists, since any valley between and b must 
lie below iria^), otherwise there would also be a peak somewhere strictly between and b 
higher than TT{ai^) which contradicts to the choice of b. Now we simply choose c to be the 
bottom of the {k, + l)-th up-stair. The proof of the decreasing statement is similar. 

Now f{a, 1) = a + 1 is proven as follows. Suppose that there are two peaks in decreasing 
position; ai < 02 and vr(ai) > 77(02), then it gives rise to a decreasing 2-staircase by the 
observation above. Hence if we have a + 1 peaks, and there is no decreasing 2-staircase, then 
the peaks must be monotonically increasing, which gives rise to an increasing (a + l)-staircase 
inductively applying the observation above. The relation = (3 + 1 is analogous. 

The proof of ( |B.1| ) is done by induction on a + f3. Consider a permutation that has at 
least f{<y,(3 + 1) peaks. We have to show that it either has an increasing {a + l)-staircase or 
a decreasing (/3 + 2)-staircase. 

The permutation up to the first f{a — + 1) peaks either has a decreasing (/? -|- 2)- 
staircase or it has an increasing a-staircase. In the latter case let ai < 02 < . . . < 0^ be 
the location of the tips. Now we look at the next f{a,(3) peaks. If they have an increasing 
(a + l)-staircase (within themselves), we are done, so can assume that it has a decreasing 
(/? -|- l)-staircase, let 61 < 62 < • • • < be the locations of the tips and certainly aa < bi. 
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Now if 7r(aa) < 7r(6i), then we found a new peak after and above 7r(aa), hence by the 
observation on increasing peaks there is an increasing [a + l)-staircase. 

If 7r(aa) > 7r(6i), then we found a higher peak before 7r(6i), hence by the observation on the 
decreasing peaks, there is a decreasing (/3 + 2)-staircase. This shows that if the permutation 
has at least f{a — 1, /3 + 1) + /(a, P) peaks, then it either has a (/3 + 2)-staircase or it has an 
increasing a-chain. □ 
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